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FOREWORD 


The objective of Contract NAS8-32408 was to assess the characteris- 
tics required by the Data Systems of future OSTA Earth-Observation systems 
of the next decade. 

It specified investigations in the six application areas Agriculture, 
Water Resources, Land Use, Global Weather, Severe Storms and Marine 
Weather. Reported here are the results of the analysis of Agricultural 
Data Systems. The authors wish to express particular thanks to Mr. Harvey 
Golden and Dr. George McDonough of Marshall Space Flight Center for their 
enlightened leadership: and to Dr. Charles T. Faludan, Dr. Sanford W. 

Downs, Jr. and Dr. Robert Jethro of Marshall for their many valuable sug- 
gestions. 

Significant portions of this work were developed in prior efforts 
under the leadership of Dr. Louis Walter and Mr. William Stroud of 
Goddard Space Flight Center, and Mr. Leonard Jaffe of NASA Headquarters. 
The author's thanks extend to them also. 

Dr. James C. Harlan provided unstinting cooperation in the difficult 
task of analysis of the spectral data bank gathered by his organization 
in Bushland, Texas. 
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1.0 SUMMARY. RESULTS AND CONCLUSIONS 


The theme for the Data System task In Agriculture, chosen In concert 
with the Marshall Space Flight Center, was "Global Crop Production Fore- 
casting." Selection of a specific theme within the broad area of Agri- 
culture was motivated by the desire to lend focus to the task; simul- 
taneously, Marshall's directives were to provide results sufficiently 
general as to be applicable to other agricultural uses of remotely sensed 
data. 


The analysis herein reported applies to a generic agricultural ob- 
servation system: LAMDSAT D was used as a benchmark against which to 

compare quantitative results. 

Two principal scenarios motivate Global Crop Forecasting: "Promotion 

of the US Balance of Trade" and "Promotion of Global Freedom from Want." 
The former was chosen as the more realistic motivator of US efforts. 

The crops to be surveyed were the major US export staples. They 
were selected from the cumulative ten-year export history of the US. 

The rationale for selecting countries candidate for survey excluded non- 
producers of specific target crops: It Included overproducers , i.e. com- 

petitors to the US for exports, and marginal producers which steadily or 
periodically import US agriproducts. 

The selection methodology yielded the staple crops wheat, soybeans, 
com, cotton, tobacco, rice, grain sorghum, and 24 countries growing 
these commodities. Of these crops, tobacco and cotton were excluded 
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from further consideration because they do not constitute alimentary 
staples. 

Substitution of crops was found to be negligible in human food, 
Important in animal fodder. This added to the list of commodities the 
small grain fodder crops barley, oats and rye. 

Crop production is the product of the two essentially Independent 
factors of acreage and yield. Since measurement of the latter from re- 
motely sensed data is as yet in the experimental stage, the study con- 
centrated on the measurement of acreage. 

The user requirements were assessed as being approximately equal to 
the accuracy performance of current conventional crop acreage Inventory 
systems. The corresponding systems specifications for accuracy range 
from 98% for US wheat to 92%, one sigma, for foreign forage crops. 

The performance of OSTA data systems in acreage mensuration was as- 
sessed from analysis of the performance achieved by ERTS and LANDSAT ex- 
perimenters. Average Inventory accuracy was found to be of order 80% in 
the time frame 1972-74; Improving to approximately 86% between 1974 and 
1978. These levels of performance were found to fall short of the user 
requirements by a factor of at least 2, at most 5. 

This finding led to investigate the data system from end-to-end, in 
order to seek the reasons for the shortcoming, and to determine which 
data system drivers required improvement. 

The analysis began at the data systems front'Cnd. i.e. at the 
interface between the natural phenomena determining crop spectral re- 
flectance and the sensor. To this effect, several of the fine-grained 


calibrated spectra from the LACIE supersltes were analyzed in detail* 

It was concluded that unique crop spectral signatures do not appear to 
exist; that the values of the reflectances fluctuate significantly * even 
within the same field; that the information contained in the differential 
amplitudes of spectral bands is in and by itself insufficient to achieve 
discrimination accuracies satisfying the user requirements; and that the 
driver of discrimination is the relative average amplitude, or contrast 
ratio, between the mean spectral amplitudes, brough about by the marked 
differences of mean reflectance which accompany the natural stages of 
phenologlc development. 

These findings led to the development of a theory of radiometric 
discrimination employing the mathematical framework of the theory of 
discrimination between scintillating radar targets. The theory indi- 
cates that the functions which drive accuracy of discrimination are the 
contrast ratio between targets, and the number of samples, or pixels, 
observed. 

The theoretical results led to three primary consequences as regards 
the data system: 1) agricultural targets must be imaged at correctly 

chosen times, when the relative evoluatlon of the crop's development is 
such as to maximize their contrast; 2) under these favorable conditions, 
the number of observed pixels can be significantly reduced with respect 
to wall-to-wall measurements; 3) the remotely sensed radiometric data 
must be suitably mixed with other auxiliary data, derived from external 
sources . 

Analysis of the optimal timing requirements in the presence of 
naturally occurring cloud cover led to the conclusion that the current 
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configuration of earth-observing platforms operating in the optical and 
IR spectral ranges need modification as to orbits and swath widths, and 
could probably substantially benefit from the addition of active micro- 
wave sensors. Analysis of the optimal sampling requirements indicated 
that the number of pixels imaged and processed drives the data system's 
operating costs: the use of a sampling scheme In heir of the current 

wall-to-wall imaging mode could reduce the system's recurring costs by 
up to a factor of two. Analysis of the role of the auxiliary data in- 
dicates that a moat important contribution to the system's price-perfor- 
mance would be a world-wide, uniform, photo graphic- type coverage of re- 
latively good resolution. Such coverage need not be provided In real 
time: its refresh rate could range upwards of three years. These 

findings can form the basis for considerable enhancement of the price/ 
performance of future OSTA agricultural observation systems. In parti- 
cular, the RBV could become the source for the auxiliary photocoverage. 
Additional savings are possible by on-board processing of the data, con- 
sisting of the addition of short-term storage. 

An additional systems savings can result from relaxation of the 
multlspectral scanner's sensitivity. This la brought about by the pre- 
sence fluctuations, found to exist around the reflectances of agri- 
cultural targets. Within the overall system's error budget, of which 
the "scene fluctuation" is the primary driver, sensor NE Ap's of order 
64 levels appear adequate. 

The resulting 4:1 reduction of sensor sensitivity can be traded off 
against lower-cost sensor implementation: or, it can be exploited to 

increase the sensor's geometric resolution. It can be shown 
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theoretically that the tradeoff In favor of performance would Improve 
LANDSAT D's resolution to approximately 15 meters. 

The theoretical results were simulated on the Marshall Data System 
Dynamic Simulator (DSDS). The results confirmed the need for modifying 
the deployment of earth-observing platforms* and the savings accruing 
from observing a reduced number of sample picels. The simulation neither 
^ confirmed nor disproved the theory* because the DSDS had not reached Its 

full Stage 111 configuration* thus could not* at the time* accept inputs 
from live scenes. The theory was however tested Independently by per- 
forming regional crop inventories In selected regions of the world* with 
good agreement between results and predictions. 

1 • 


2.0 SCENARIO 


2.1 Significance of Global Crop Survey and Forecast 

The activities of estimating and forecasting agricultural production 
are of major Importance virtually everywhere In the world. 

Inaccurate forecasts cause agricultural producers to make erroneous 
production decisions: they distort optimal Inventory carryovers. The net 

result is lowered social efficiency through wasted efforts. Improved fore- 
cast accuracy reduces this waste motion > hence produces a net increase In 
social welfare. 

Precise forecasts obtained sufficiently early permit the formulation of 
national plans and policies best suited to cope with surplus or scarcity. 
Good predictions allow agricultural producers to take remedial actions 
aimed at Increasing output, reducing costs, or taking shelter against de- 
clining prices caused by overproduction. The earlier this Information is 
available, the greater the possible spectrum of ameliorative actions and 
the larger the potential benefits. 

The worth of agricultural survey and forecast activities is evidenced 
by the efforts and resources which are being devoted to It worldwide. 

Table 2-1 presents a conservative estimate of the expenditures by world 
governments for agricultural data gathering. Additional moneys are ex- 
pended yearly by commodity dealers, traders. Chambers of Commerce, Inter- 
national organizations, other private and public enterprises [1]. 




TA81E e-1 


ycwa.y HORinmoc expbpituiies fm chop suRVE»/Fo«EC/gT »a 77 > 


region 

SURVEY/FORECAST COST 
MILLION U.S. S. 1977 

U.S. (OOHESTIC < FOREIGN) 

40 

CANADA 

2.4 

WESTERN EUROPE 

48.2 

EASTERN EUROPE 

8*9 

LATIN AMERICA 

2*5 

AFRICA/NIO-EAST 

2*1 

ASIA EXa. USSR, PRC 

4*5 

OCEANIA 

1*5 

WORLD 

no.i 


The world's important crops number in the hundreds. Most Important in 
terms of market value, quantity and caloric content are the staple crops. 

Figure 2-1 shows that grain crops are the most significant: they pro- 

vide over half of mankind's food. [2j , [3]. 

2.2 World Scenarios 

The configuration and characteristics of a Data System for global 
agricultural forecast are a function of the surface area, geographic loca- 
tion, number of crops to be surveyed. 

Two principal scenarios determine the selection of the countries and 
crops candidate for survey: 





Figure 2-1 


PANORAMA OF THE WORLD’S PRODUCTION 
OF STAPLE CROPS 



BY WEIGHT 



BY VALUE (AT INTERNATIONAL EXPORT PRICES) 
WCY's 



BY CALORIC CONTENT 



1 


As seen from the U.S. national outlook 


2. As seen from the outlook of the international community, l.e, 
the United Nations. 

The U.S. national goals have essentially remained unchanged throughout 
U.S. history. They are: 

f Maximize the citizen's well-being 

• Foster the U.S. balance of trade 

• Optimize U.S. foreign policy 

The United Nations goal is essentially the humanitarian objective: 

• Foster world peace 

• Promote global freedom from want 

2.3 The U.S. Scenario 

In the agricultural arena, the three U.S. goals are interrelated. 

Promotion of the balance of trade Implies essentially "sell more, at 
more advantageous prices." 

Citizen's well-being, or consumer satisfaction, means: "but without 

unduly increasing domestic prices for U.S. -produced commodities, and main- 
taining favorable prices for Imported agricommodities." 

The opportunity for a significant foreign policy role arises if the 
U.S. can trade food in exchange for political or socieconomic considerations. 

This study focused on the first goal, l.e. maintenance and promotion 
of a favorable balance of trade. 
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2.4 Derivation of the Crops to be Surveyed 

Agricultural products have always played a major role In the U.S. 
export picture. In the mid-nineteenth century they represented approxima- 
tely 80% of total U.S. exports. Although the agricultural share of total 
U.S. exports has since declined to 25%, the gross dollar value of agricul- 
tural exports has Increased from $150 million In the 1870's to over 
$25 billion, In current dollars. 


U.S. agricultural exports fluctuate yearly. In response to domestic 
and foreign alternations of production shortfalls and surpluses. To smooth 
out these f lucutatlons , the export values should be Integrated over a suf- 
ficiently long period of time. Table 2-2 ranks the U.S. agriexports based 
upon their 10-year cumulative value [3]. 



TABLE 2-2 




U.S. GROSS AGRICULTURAL EXPORTS. 10.VEAR CUMULATIVE VALUE 


19U-1974, 

CURRENT $ HIUION. 

FOB 



COmODlTY 

ABSOLUTE 

VALUE 


RELATIVE 

VALUE 

WHEAT 

$10*527 


16X 


SOYBEANS (FOR FEED) 

14«750 


1SX 


CORN 

13,873 


14X 


AltlHAL A ANIMAL PROMCTS 
OF MUCH 

10*802 


11X 


FATS OILS A GREASE 


2*603 


(3) 

HEAT A HEAT PRODUCTS 


1.8GG 


(2) 

OIL SEEDS 
OF WHICH 

9*090 


9X 


SOV OILCAKE A NEAL 


4,299 


(4) 

COTTONSEED OIL 


660 


(1) 

COTTON 

6*525 


6X 


TOBACCO 

6*186 


6X 


FEED GRAINS* ExaUOiNG CORN 
OF WHICH 

4*021 


4X 


GRAIN SORGHUH 


3*004 


{3) 

RICE 

3.664 


4X 


SUBTOTAL 

87.548 


B6X 


OTHER 

14.481 


14X 


TOTAL 10-TEAR EXPORTS 

$102*029 


tool 




Table 2->3 syntheslzea the Information of Table 2-2. It excludes the 


animal products* for the following reasons: 


• The consistency of the animal Inventory is not directly amenable to 
survey from Satellite Remote Sensing. 

• The diet of animals is composed of grains* pasture and forage 
in varying proportions in differing countries. The feed grains 
are already included in the listing of Table 2-3. 


TMIE i-3 


CROSS U.S. EXroitTS OF STAflE CROPS. 1Q.YEAR CUHUIATIVE VALUE 
19M-1974. CURRENT S MILLION. FOB 



V/^IUE 

PERCENT 

sov m soy rrooucts 


21 

MEAT 

18,527 

20 

CORN mo CORN PRODUCTS 

14,870 

16 

COTTON 

6,525 

7 

TOBACCO 

6,185 

7 

RICE 

3,684 

4 

SORGHUM 

3.004 

3 

SUBTOTAL 

71,845 

79 

OTHER LESS MVWL 
PRODUCTS 

19.292 

21 

TOTAL LESS ANIMAL 
PRODUCTS 

91.137 

100 


Percents tony not a4d due to round1n9 errors. 


Rangeland (pasture) and forage (hay) are also important candidates 
for survey for two reasons: 1) they are complementary or substitute 

animal foods* and can thus affect the demand for the primary grains; and 
2) they act as "confuser" crops in the remote observation of the primary 
Staples . 
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Table 2-3 indicates that soybeans, wheat and corn head the list of 
exports by value. They constitute 58% of U.S. staple agriproducts exported 
during the decade, exclusive of animal products. The other major categories 
combined - cotton, tobacco, rice, sorghum - represent 21%. The aggregate 
of the several minor commodities is 21%. No single commodity below the ranking 
of sorghum produced exports greater than 1% of the total. 

Before concluding that the seven crops indicated in Table 2-3 are 
indeed those of major interest to U.S. agricultural exports, it is well 
to analyze whether the 10-year cumulation does not smooth out and hide 
major fluctuations. Figure 2-2 Indicates that this is not the case: 
the relative percentages of the exports of the staple crops vis-a-vis 
the total agricultural exports have remained relatively constant over the 
last decade. Within the validity of historical extrapolations, and barring 
major exogenous disturbances, they can thus be assumed to remain substan- 
tially the major U.S. agriexports of the near future. 

2.5 Derivation of the Countries to be Surveyed 

2.5.1 Logic 

The "universe" for export of U.S. agricultural commodities subdivides 
as follows: 

1. Coiuitrles which do not consume the product at all: i.e., non- 

customers. 

2. Countries which are consistently self-sufficient (for a given 

agriproduct which the U.S. produces in surplus): neither cus- 

tomers nor competitors. The sub— category of Countries 
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which do nofc purchase from the U.S, for policy reasons, valid 
in the 1960's, has, since the massive entry-to-market of USSR, 

PRC and Eastern Europe, dwindled into insignificance, at least 
for the time being. 

3. Countries which produce the product, but not in sufficient* amounts 
to meet their demand, and which are steady importers. 

4. Countries wherein the difference between production and internal 
demand alternates between surplus and deficit. These are poten- 
tial U.S. customers in shortfall years, potential competitors In 
surplus years. 

5. Countries which are consistent exporters. These are potential 
competitors of the U.S. 

6. The U.S. Itself, to assess its year-to-year capability to 
export . 


The requirements for agricultural surveys by the U.S. vary among these 
categories of countries. 

Countries which exhibit characteristics 1 and 2 - l.e. consistent 
non-customers and non-competitors for any given product - are not Important 
to agricultural U.S, exports of that product. Surveys of their agricul- 
tural production of that product are thus of little interest to the U.S. 
foreign trade posture. 

Surveys of countries in category 3, i.e. , which exhibit the qualities 
of marginal producers and steady and consistent Importers are of Interest 
to assess the quantity of product they are likely to Import in any one year. 
This assessment serves as the basis for U.S. agriexport planning on how to 


convert this "addressable market" into sales for the U.S. The interest is 
proportional to each country's volume of the addressable market and to the 
potential competition from other surplus-producing Countries* 

Countries of category 4» fluctuating customer-competitors » are of 
Interest in proportion to the magnitude of the fluctuations of their ag- 
ricultural production. 

Surveys of countries of Category 5« the potential competitors, are of 
interest in proportion to the magnitude of their surplus. 

Figure 2-3 depicts the logical choices. 

2.5.2 Procedure for derivation of the countries to be surveyed 

The procedure Is to isolate, down to a meaningful level of significance, 
the countries which fall In categories 3, 4, and 5, for the principal crops 
of Interest to U.S. agriexports. 

Figure 2-4 depicts the principal customers for the major U.S. export 
crops over the last decade. They are ranked from bottom to top by percen- 
tage of the dollar volume of Imports from the U.S.; from left to right by 
total dollar value of each imported commodity. The blank unlabeled areas 
at the top of the histogram bars comprise a rather large number of customers, 
each of which Is significant at most to the 1 % level; l.e.. Its Imports 
during the decade did not exceed more than 1 % of the total U.S. exports of 
the corrspondlng commodity. 

Note that Figure 2-4 Indicates percentages and not absolute dollar 
values. The dollar purchases of each customer can be obtained by 
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multiplying the percentages of product acquired* from Figure 2-4* by 
the total dollar value of the product exports* from Table 2-3. For 
example* Japan's purchases of wheat were 10% of the total U.S. exports; 
those of sorghum* 45%. The total U.S. exports were $18,5 Billion for 
wheat* $3 Billion for sorghum. Thus Japan's 10% share of U.S. wheat 
exports brought a revenue of $1.85 Billion: Its purchases of sorghum* 

in spite of the much higher share* resulted In a lesser revenue* e.g. 
$1.35 Billion. 

The Information of Figure 2-4 allows elimination of the countries 
which fall In categories 1* non-consumers* and 2* self-sufficient. It 
aggregates categories 3* net Importers* and 4* fluctuating importers- 
exporters. 

The net importers fall Into two sub-categories t steady Importers* 
typified by Japan, whose purchases essentially repeat year after year; 
and fluctuating Importers. This second sub-category is clearly of 
greater Interest to U.S. survey efforts. Figure 2-5 exemplifies Its 
behavior. 

Table 2-4 serves to assess which Importing countries fall into 
category 4* fluctuating customer-competitors. It Indicates that the 
USSR Is the only such for the commodity wheat, since during the decade 
It did export as well as Import. The PRC likewise Is a fluctuating 
customer-competitor for the commodity' soybeans* France for com. 
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FICURE 2-5 

THE MOST FLUCTUATING U,5. CUSTOMERS 
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TABie 2-4 

WJOR HOBLD tMPOBTS > EXPORTS OF StlECTEO AGRlCtHTUMl. COHHDOiTlES 
6»YEAf> CUHUIATIVE VALUE (IWIB-yBl 
«H CUAREHT - gpORTS FOB, IMPORTS CIF 

WHEAT com SOtBEAHS 


WORLD 

EXPORTS 

26.830 

IMPORTS 

28.1S0 

WORLD 

U.S. 

9.740 

- 

U.S. 

CANADA 

6,050 

* 

FRANCE 

ffUMCC 

3.160 

- 

ARGENTINA 

USSR 

2.730 

2.040 

SOUTH AFRICA 

AUSTRALIA 

2,150 

* 

JAPAN 

ARGENTINA 

840 

- 

ITALY 

PRC 


2.430 

UK 

JAPAN 


:70 

FRG 

UK 


2,060 

NnHERLANOS 

INDONESIA 


1.500 

SPAIN 

FRG 


1,400 

USSR 

BRAZIL 


1,040 


ITALY 


840 


EGYPT 


740 


CZECHOSLOVAKIA 


610 


PAKISTAN 


310 



EXPORTS 

INPORTS 


EXPORTS 

IMPORTS 

13,600 


WORLD 

9.620 

9.720 

7,no 

- 

U.S. 

0,440 

- 

1,640 

150 

BRAZIL 

706 

w 

1,520 

• 

PRC 

350 

500 

660 

* 

JAPAN 


2.590 


2.540 

FRG 


1.620 


1,990 

NETHERLANDS 


S30 


1,380 

SPAIN 


910 


1,370 

ITALY 


660 


1,330 





970 





630 





Table 2'-4 also Identifies the countries falling within category 5, the 
steady competitors. These are countries which do not import because they 
possess a consistent, repeating discretionary surplus; i,e, their pro- 
duction consistently exceeds Internal demand. 

Figure 2-6 presents a panorama of the principal competitors. It 
indicates that the U.S. not only possesses the largest but also the 
most consistent and reliable surplus. 


2.5.3 Recapitulation of Crops and Countries Candidates for Survey 

Table 2-3 recapitulates the principal crops and countries candi- 
dates for survey of the major staples for human consumption, down to 
the 2% level (as percentage of total U.S. agriexports) for the crops. 


Figure 2-6 

DISCRETIONARY SURPLUSES OF MAJOR AGRICULTURAL EXPORTERS 
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and to the 2% level (fraction of imports of a given commodity from 
the U>S.) for the Countries. Tobacco is not included since it is not 
a major staple for human consumption: it ls» rather, a luxury Item. 


T/miE 2-5 


CROPS AHO COUNTRIES IKITIAL CAHDIOATES FOR SURVEY 


WHEAT 

CORN 

SOYBEANS 

RICE 

COTTON 


U,S.5*R* 

- 

- 

U.S.S.R. 

U,S, 

U,5. 


u.s« 

U.S, 

PRC 

PRC 

PRC 

PRC 

PRC 

INDIA 

INDIA 

INDIA 

INDIA 

INDIA 

FRANCE 

FRANCE 

- 

• 

- 

CANADA 

- 

- 

“ 


AUSTRALIA 

- 


- 

“ 

TURKEY 

* 


- 

TURKEY 

ITALY 

ITALY 


ITALY 

- 

POLAND 

- 


- 

- 


BRAZIL 

BRAZIL 

BRAZIL 

BRAZIL 

. 

S. AFRICA 


BRAZIL 

S. AFRICA 

ROMANIA 

ROMANIA 


- 

- 

- 

- 


INDONESIA 


FRG 

. 


- 

- 

YUGOSLAVIA 

YUGOSLAVIA 


- 

• 

ARGENTINA 

ARG^TINA 


- 

ARGENTINA 

HQROCCO 

- 


- 

• 


HUNGARY 


- 

- 


EGYPT 


- 

EGYPT 

- 

- 


bakgladesk, 

JAPAN* THAILAND 

" 

. 



* 

PAKISTAN 


_ 


- 

MEXICO 


Before concluding that this is Indeed the final set of crops and 
countries, it is Important to Investigate whether other crops produced by 
foreign customers are used as substitutes for the primary V.S. exports: 
for, in this case, and to the extent of the possible substitution, foreign 
crops other than those indicated in Table 2-5 may become candidates for 




s urvey . 
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2.5.4 The Effect of Substitution In Hmnan Nutrition 

Past studies and field surveys by U.S. Agencies and Foundations, by 
the UN and the FAO, have uniformly concluded that substitution, even 
of very similar foods, is not a viable alternative for human consumption 
even in periods of famine. At least, not In the short-term; habits 
and tastes change but slowly £43. 

This is confirmed by Figure 2-7, drawn from FAO statistics, which 
depicts the historical evolution of human consumption of cereals in 
selected countries in the three primary categories of wheat, eaters, 
maize eaters, and rice eaters. Because the per-caplta income has grown 
steadily and significantly over the years shown, the curves illustrate 
the well-known relationship that ttie consumption of staples grows 
initially with increasing per capita Income: beyond a threshold level, 

further Increases in income result in a decreasing pro-capita consumption 
of staples in favor of protein and "Itixury** foods. More significant 
to the question of substitution. Figure 2-7 indicates that the ratios of 
the diverse staples consumed by a given country change but slowly. Also, 
short term reductions in pro-capita Intake of one staple are not 
compensated by corresponding Increases in other staples. Ihls indicates 
that siibstltution is not an important driver in human nutrition, except 
over relatively long time spans . 
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2 . 5,5 The Effect of Substitution In Anltnal Feed 


Over the last fifteen years« the Increased demand for protein has 
caused a significant shift in the composition of animal feed. Appendix A, 
Figure A-4 Illustrates the trend for the USSR and the U.S.; whereas human 
grain consumption has remained substantially constant* that of feed has 
significantly increased. 

Currently* of the world's cereal harvest only approximately 37% Is 
dedicated to human consumption* the balance 50% after deducting allocations 
for seed* waste* Industrial usage. Is devoted to animal consumption. 


Contrary to what happens In htiman consumption* cereals of all types 
can be substituted in animal feed. As such, their survey Is Important 
within the customer and customer-competitor countries already selected 
In Table 2-3* and in other countries which are major producers of these 
commodities. The ranking of the major producers of feed grains is pre- 
sented in Figure 2-8. Cotton is also included. 


The complete tabulation of crops and countries which are the final 
candidates for survey is presented in Table 2-6. 

It is well to reiterate that these crops and countries are only those 
that have passed the criteria established previously. Japan, for example, 
although a major and steady customer for several US agricultural staples, 
is not a significant producer of these: as such* it Is excluded from the 

list of candidates for survey* simply because there would not be much to 
survey. 

1 -; 
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FIGURE 2-8 

RANKING OF WORLD PRODUCERS OF PRINCIPAL CEREAL SUBSTITUTION AND FIBER CROPS 

SOURCE: USDA AGRICULTURAL 
STATISTICS 
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TABLE 2-6 


WHEAT 

U.S-S.R. 

U.S. 

PRC 

INDIA 

FRANCE 

CANADA 

AUSTRALIA 

TURKEY 

ITALY 

POLAND 


ROMANIA 


FRG 

Yt^SLAVlA 

ARGENTINA 

MOROCCO 


CROPS AND COUNTRIES Fm CANDIDATES FOR SURVEY 
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2.6 Effect of Global Crop Forecast on Other U.S. Goals 

2.6.1 Effect of the requirement to maximize consumer satisfaction 

In a free-market economy, early prediction of the status of crops 
will restilt in the maximum possible reduction of U.S. consumer prices 
for the primary commodities. The Ut.)S> spends approximate one-half of 
its revenues from the exports of agriproducts in importing complementary 
and suppementary agricultural commodities and derived products. 

Complementary Imports are items which are not produced in the U.S,: e.g. 
coffee, bananas, tea. Supplementary Imports consist of products 
available domestically, but which consuners cherish for reasons of taste or 
custom: e.g. wines, canned tomatoes, pasta products. 

The complementary and supplementary agricultural products play a 
significant role in insuring the satisfaction of U.S. consumers. Thus 
inventories of the respective producing countries - for example, 

Brazilian coffee - are important, and a scenario analogous to, but in 
reverse of that presented above for the exports, could be constructed. 

This was, however, beyond the scope of this study. 

2.6.2 Role of Global Crop Forecast as Tool of Foreign Policy 

The optimum trading of U.S. -produced food in exchange for political 
advantages requires the early knowledge of the requirements of the numeroiis 
countries who are actual or potential AID recipients. This implies that 
the survey of principal crops should be extended to at least the most im- 
portant of these countries. The rationale for selecting target crops 
and countries differs from that developed above: it addresses the selection 

objectives of U.S. Foreign Policy. Crops different from those thus far 
identified do emerge — e.g. cassava. 
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Thls investigation was beyond the scope of this study. 

2.7 The International Scenario 

The contribution achievable to the UN goal of "world peace" by 
providing adequate nutrition to the world's masses is a matter of philo- 
sophical speculation. History provides abundant examples of wars and of 
local disturbances Induced by want; as well as examples of prolonged mass 
deprivations not leading to such violent epilogues; and also Instances of 
major wars initiated by well-fed contenders. 

The second UN goal, "promoting global freedom from famine" is more 
specific and amenable to quantitative analysis, without the need to pass 
judgements on its philosophical correctness. It is presented in 
Appendix A. 

The UN/FAO scenario was not utilized for this study. It Is however 
worth noting, as shown in Appendix A, that this scenario maintains the 
crops and countries of the U.S. Balance of Trade scenario, except for 
the Incltislon of additional countries producers and consumers of the 
staple crop rlre. 
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3.0 SYSTEM PERFOBMANCE REQUIREMENTS AND CURRENT PERFORMANCE 


3.1 Derivation of the Data System's, Performance Requirements 

Appendix B presents an overview of the agricultural survey aethods 
currently In use throughout the world. As Indicated therein, crop 
production forecasts comprise two distinct functions: 

• Measurement of the current status of the crop. For reasons of 
cost, this measurement Is always performed by statistical sampling. 
Following the USDA designation. It Is also called "estimate." 

• Extrapolation to the future harvest, more properly termed "fore- 
cast." This function Is accomplished based nn historical regres- 
sions augmented by "best estimates" of the future Incidence of 
adverse or beneficial phenomena affecting the crop, such as 
weather. 

The estimate of the current status Is composed of two subfunctions; 

• Estimate of the acreage under cultivation 

• Estimate of the crop's condition. 

This study focuses upon the acreage estimates, because they are 
currently amenable to measurement from orbital data systems. 

Table 3-1 presents the accuracies of acreage mensuration achieved by 
various countries and by U.S. States and Counties [5], 16], [7]. It Is 
subject to the following Interpretation: 
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• The reported accuracies were obtained from direct discussions with 
responsible individuals in the US and abroad, augmented by analysis 
of their data. 


• The reported accuracies are valid late in the season; i.e. not 
earlier than at least one month prior to harvest. 


• The stated accuracies refer to sampling errors only, i.e. the com- 
bination of the random errors Induced by the limited number of 
samples and by the accuracies of measuring each sample's crop 
acreage. Systematic errors are not included. 



TABLE 3-1 


ACCURACIES CWAENTLV ACHIEVED IH THE MEASUBEHEWT OF ACREAGE 


COUNTRY 

CROP 

ACREAGE HEASUREMENT 
ERROR PREHARVEST 
ONE SIGMA 

o.s. 

WHEAT 

12: 


CORN 

16 : 


SOY 

16: 

ITALY 

WHEAT 

14: 


FORAGE 

ib: 


CORN 

i71 

WESTERN EUROPE 

SIMILAR TO ITALY 


CANADA 

SPAING WHEAT 

13X 

USDA-FAS WORLD 
ESTIHATE 

WHEAT 


0*S* STATES 
(without auxi- 
liary list 
sampling) 

HAJOA STATE CROP 

-V17S-15: 

U.S. COUNTIES 
(without auxi- 
liary Ifst 
sampl Ing 

MAJOR COUNTY CROP 

UP TO 120-25S 




i 
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Perusal of Td>le 3-1 leads naturally to the question as to whether 
the accuracy figures therein presented do Indeed represent the user 
requlreiuentB. In numerous discussion by EOOsys terns personnel with do- 
mestic and foreign user agencies, the following transpired: 

• Currently achieved accuracies obtained at tkt c.uM.ent itatutooy 
Aapo^vting datu cannot be reduced. 

s They must be achieved at costs not in excess of current costs: 
in fact, it is desirable that the costs be lower. 

s Accuracies lower than those currently achieved are only useful 
if they can provide valuable adjuncts to the current conventional 
systems, either in terms of significantly lowering costs or of 
providing significantly Improved performance of the total system. 

a Users are willing to pay a modest premium for significantly improved 
accuracies achieved at the. &tatxxtofiy fLepoftiing datz&, 

m Users are disposed to compromise accuracy for timeliness, l.e. 
less accurate results obtained earlier than currently achieved 
are acceptable. 

S The combination of accuracy and timeliness warrants a significantly 
higher premium than what users are willing to pay to accuracy alone. 

The current "user market” sitiiatlon strongly counsels that the ac- 
curacy levels depicted in Table 3-1 be assumed as goals for the data 
system's performance. 
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3 . 2 AssesBinent of the Performance of Orbital Data Systems In Estimating 

Ac r sage 

The purpose of this Investigation was to quantify numerically the 
performance of existing orbital Data Systems, for comparison against the 
user requirements depicted In Table 3-1. The Investigation was performed 
In two steps: 

• First, analysis of all available reports on cr<.^p acreage mensu- 
ration by ERTS Investigators; these represent tests performed 
during a concentrated activity In the early phases of the LANDSAT 
program, from its launch In July of 1972 until 1974. 

• Second, analysis of findings from later tests, to assess the 
trend of improvement. 

To this effect, results were gathered from available, published ERTS 
Investigation reports containing niunerlcal performance figures. Details 
of the analysis are presented in Appendix C. 

The resulting statistical population consisted of 224 distinct tests , 
performed in 10 locations, embracing 7 crops. 

The cumulative curve of results Is presented in Figure 3-lA. It 
graphs the accuracies achieved versus the percentage of the tests achieving 
them. It presents the results for both the inventory and the mapping, or 
land use, mode. Inventory mode means the determination of the proportion 
of a given crop within the area under examination, without however speci- 
fying the crop's geographic location. It Is the mode universally used for 
estimating agricultural acreage. Mapping mode, as the word Indicates, 
connotes the identification not only of the species sought, but also of 
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FIGURE 3-1 

DISTRIBUTION OF CUSSIFICATION ACCURACIES FOR ERTS INVESTIGATIONS ON 

ACREAGE MENSURATION 
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thelr geometric boundaries and geographic locations > This mode is gene- 
rally not employed in conventional agriculture^ acreage estimation, due 
to its higher cost with respect to the inventory mode. It is presented 
in this report because of its implications to the theory of discrimina- 
tion discussed in Chapter 5. 

It can be noted from Figure 3-1 that the accuracies achieved by ERTS 
investigators in the mapping mode are lower than those competing to the in- 
ventory mode: this is quantitatively consistent with the fact that the 

mapping mode requires more Information than proportion estimation: con- 

versely, for equal information content, its level of performance should 
be lower. 

Since, as will be shown in Chapter 5, the accuracy of multlspectral 
classification worsens when the proportion of target crops to the total 
scene is small, a subset was extracted from the total statistical popula- 
tion, based on the criterion of 10% or larger proportion of target crops. 

This yielded 64 distinct data points: the test results are presented 

in Figure 3-lB. Note the higher accuracies achieved. 

The probability density functions corresponding to the curves of 
Figure 3-1 are shown in Appendix C. Computation of the accuracies as a 
function of confidence limits from these yields the classification re- 
sults depicted in Table 3-2. It shows that the one sigma accuracy of 
classification in the in\>Qi\toKy modo. for the total population of 224 
data points is 74%: this Increases ' to 86% for the subset of 64 data 

points featuring 10% or greater proportions. 
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ACCURAcy OF CLASSIFlOTtOH OF CBOPS *CHltVED 
»y EATS INVEST 1CATOR5 



Uftivftrse of 224 cUssifIcatlon tests on 7 crops 
one confidence 

• 10 geogrephlc locetlons. 


PROPORTKM 

ESTtHATtON 

HAPPir^G 

ACCURACY : 

74 

63 

imn 1 

26 

37 

64 cU&sU1c«t1qn tests extracted front the above uni verse » wherein the 
crop‘to-tota1*«rea proportion greeter then 10S. 

ACCURACY X 

66 

76 

CtlROR % 

14 

24 


Comparison with Table 3-2 shows that the average performance of ERTS 
investigator tests in accomplishing Inventory fell short of the accuracies 
currently achieved by conventional methods, with the possible exception of 
inventories performed at county levels, 

ERTS investigations ceased after 1974. To assess whether the early 
performance has improved in the approximately four years elapsed since, 
a number of subsequent investigations, performed by reputable organiza- 
tions, were analyzed. The details are presented in Appendix C; the 
overall results are shown in Table 3-3. It can be seen that there has 
been improvement since 1974; yet performance is still short of the 
user's requirements. 


-36- 




TABLE 3-3 

ACCWACY OF CLASSIFICATIOW OF CROPS ACHIEVED FROM 1.974 TO 1978 

UnWcrse of 16 cUssITIcatlon endiavon on 10 craps, 13 googriphlc 
locations; one alffna confidence. 


ACCUMCT X 
ERROR X 


PROPORTION 

ESTIHATIOH WPP1H6 


66 

73 

}A 

27 


The logical next inquiry is to ascertain wheth'»r the 75-85% inven- 
tory accuracy thus far achieved by LANDSAT experimenters is a limit im- 
posed by nature, or whether there exist driving elements in the Data 
System whose exploitation can further enhance crop discrimination and 


Identification . 
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4.0 INVESTIGATION INTO THE PHYSICAL REASONS FOR THE DISCREPANCY 
BETWEEN REQUIREMENTS AND PERFORMANCE 

4.1 Methodology 

The generalized end-to-end data system of radiometric earth-observation 
platforms employed currently and expected to be used In the 1980-1990 decade 
Is depicted In Figure 4-1. The methodical Investigation of Its performance 
was begun at the frort-end, by analyzing the relationships between the In- 
formation provided by nature, l.e. the crop -spectra, and the measurement 
capabilities of radiometric sensing. 

The Initial question addressed was whether reflectance spectra of 
crops do provide sufficient Information for accurate Identification or 
discrimination, and, If so, under what operating conditions. By discri- 
mination is meant the recognition that the observed classes are distinct: 
identification, In addition. Is the attribution of a specific label, e.g. 
agricultural crop, to each distinct class. The difference between dis- 
crimination and Identification Is conceptually akin to that between In- 
ventory and mapping. 

Proper performance of this Investigation suggests the removal. In- 
sofar as possible, of all effects which perturb the spectral reflectance 
patterns: primarily atmospheric absorption, path radiance. Instrumental 

drifts, effects of varying orientations of Illumination and look angles. 

This requires the availability of "pure" spectra, l.e. of fine-grained 
crop spectral measurements taken under carefully calibrated conditions. 


FIGURE 4-1 
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Several spectral data bases were available in the literature prior 
to 1975. Table 4-1 lists the principal ones. 


TABLE 4-1 

PRINCIPAL SPECTRAL DATA BASES AVAILABLE IH THE UTERAHiRE 
SPONSOR PUBLICATION 


DESIGNATION 

AGENCY 

DATE 

LOCATION 

KTF 

NASA/WTF 

12/73-5/74 

HTF > BAT ST. LOUIS, HS 

LARS 

NASA 

7/68 

LARS - WEST LAFAYETTE. : 

ERIH 

U.S. AIR FORCE 

1972 

NTIS ■ SPRINGFiaO, VA 

EROL 

U.S. ARMY 

1961-1962 

NTIS . SPRINGFIELD. VA 

NESLACO 


10/66-9/67 

USDA ARS - WESLACO, TX 

KRINOV 

USSR 

1936-1938 

NRC . OHAWA, CANADA 

IITRI* 

NASA/U5AF 

VARIOUS 

NTIS - SPRINGFiaO, VA 


* 

Sunmary of data from other sources. 


Unfortunately, most of these spectra do not correspond to ideal condi- 
tions, because: 1) they are either spectra of plant leaves taken in the 

laboratory, thus not of the plants as a whole; 2) calibrations surrounding 
the taking of whole-plant spectra were not consistent. Nevertheless, 
these spectral data bases were useful as secondary, qualitative sources of 
information. 


The first set of "pure" spectra became available from Bushland, Texas 
in 1975. They were supplied to ECOsystems through the courtesy of Dr. 

James C. Harlan of the Spectral Characteristics Laboratory, Remote Sensing 
Center, Texas A&M University. Dr. Harlan was also most generous with minu- 
tious Information as to the experimental procedure followed, ground truth 





-40- 

ntapsf files and aerial Imagery. This allowed certification of the quality 
of these spectra and permitted the compensation of spurious effects which 
had crept into the data gathering procedure in spite of the experimenter's 
diligence. 


Subsequently* other data sets, also collected under calibrated condi- 
tions at low altitude* became available. The set of data used In this 


study Is Indicated In Table 4-2. 



TABLE 4-2 



CHARACTERISTICS OF THE "PURE' 

SPECTRAL DATA USED IN THE ANALYSIS 


BUSHLANO. TX 

HAND CO. f SO 

IMPERIAL VALLEY. CA 

TYPES OF CROPS 

N, R. P. (C,S) 

H, A. P, 0, B (C) 

(H, S, SB, AS. A. HE) 

NUNBER OF FIELDS 

100 

:i5 

;ioo 

SPECTRAL RAKGE, HlCm% 

OA to 2*4 

0.4 to 2,4 

0.4 to I.O 

SPECTRAL RESOLUTION. MICRONS 

0.02 from 0*4 to 1*) 

0.02 from 0.4 to 1*1 

0.0014 

NUMBER OF SPECTRA GATHERED 
PER FIELD 

;20 

:20 

:15 

CALIBRATION PROCEDURE 

OVERFLIGHT OF REFERENCE 

PANEL 9 20' INTERVALS 

REFERENCE SOURCE BEFORE A 
AFTER FLIGHT 

INSTRUMENT 

S-19IH CONICAL WHEEL SPECTROMETER 

COLLINS SPECTROMETER 

PLATFORM 

HELICOPTER 

HELICOPTER 

AIRCRAFT 

OATES OF GATHERING 

I974-7S CROP SEASON 

1976*77 CROP SEASON 

1975 GROWING SEASON 

GROUND RESOLUTION, METERS 

23.8 

23*8 

18 X 18 

aiGNT ALTITUDE, METERS 

60.9 

60.9 

610. 

H • WHEAT R - RYE P - PASTURE C - CORN S • SORGHUM 

A • ALFALFA 0 - OATS S ■ BARLEY SB ■ SUGAR BEETS 

AS • ASPARAGUS HE • HaONS 

( ) • SPECTRA AVAILABLE DURING PARTIAL GROWTH CYCLE ONLY. 



These data sets represent a major improvement with respect to the 
earlier spectra available In the literature. Nevertheless, they still 
present significant uncertainties. For example: none of the data sets 

embraces a temporal span sufficiently wide so as to contain all the 
phenologlc stages of com or sorghum. The data-gatherer's primary 



interest appears to have been the measure of wheat. In several cases, 
obvious sensor and/or processing errors have crept in. Frequently, the 
percent of crop ground cover is not unambiguously stated: thus several 

crop spectra are mixed with soil spectra in such fashion that it is dif- 
ficult to unravel the true reflectance of either. In many cases, the 
nomenclature does not correspond with the actual ground cover: fields 

planted to a certain crop which never does emerge are labeled as con- 
taining that crop. No attempt is apparent to reconcile widely disparate 
reflectances pertaining to the same field, as for example the reflectances 
of soil, labeled "bare," before emergence and after harvest. 

This state of affairs, particularly evident in the Hand County 
spectra, reqxiires considerable Judgement, and inter-comparison of spectra 
from all the data sets, in order to av^''ld the conclusion that crop reflec- 
tance spectra are essentially random. 

4.2 Results from the Analysis of the Field Spectra 

The details of the analysis of the available data sets are presented 
In Appendix D. The salient conclusions are: 


1. The shapes of the crop reflectance spectra are significantly 
similar at similar growth stages - i.e. after the crops are 
sufficiently developed to cover the underlying soil. As exem- 
plified In Figure 4-2, there is little evidence of sharply 
unique spectral structures among the various crops. 

2. The spectral shapes remain significantly similar throughout the 
crop's growth stages, provided that the crop's ground cover Is 
sufficiently dense so as to minimize interference from soil re- 
flection. 


RADIANCE (nw. cm . sr -".nm 


m 






FIGURE 4-2 


\ , 
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3. The spectra exhibit significant differences of mean amplitudes 
as a function of stage of growth. Typlcaly progressions are 
exemplified in Figure 4-3. 

4. The mean values of the spectral amplitudes exhibit fluctuations. 
This is exemplified in Figure 4-4. Fluctuations are present 
within single fields and in ensembles of fields • and occur in 
all regions of the spectrum. 

5. The fluctuations among ensembles of fields are generally larger 
than those pertaining to single fields. A possible expl£ination 
of the Intrafleld fluctuations Is the presence of discontinui- 
ties — patches* bare spots* stressed areas — within the fields; 
and the clutter effect originating from the varying orientations 
and sizes of the reflecting components of the canopy. A logically 
appealing explanation of the interfield fluctuations is that 
cropped fields do not reach the same development stage simul- 
taneously. Figure 4-5 depicts for various crops Che fraction* 

as a function of time, of the cropped surface which reaches a 
given stage of phenologlc development. It can be seen that as 
many as 60 days are required to progress from 0% to 100% of the 
covered area. Thus an ensemble of fields generally contains in- 
dividual fields at different phenologlc development* and hence 
different reflectances. The presence of stress causes similar 
effects. The available data* cf. Appendix D* suggest 0.1 as a 
typical lower value and 0.3 as a typical upper bovmd of the coef- 
ficients of variation for an ensemble of fields. 

6. Identification was not found possible from spectral signature 
or from differences in amplitude alone. It requires additional 
external information* e.g. knowledge of the phenologlc progres- 
sions prevailing within the region being observed. 


K 


fCE FOR THE CROPS CONTAINED 
>ATA SET 

10% GROUND COVER 





5 - 




PERCENT OF CROPPED AREA REACHING INDICATED STAGE 
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FICORE 4-5 

AREAL DEVELOPMENT OF CROPS AS A FUNCTION OF TIME SOUTHWEST CROP 
REPORTING DISTRICT. KANSAS. 1973 SEASON 


MAY I JUNE I JULY | AUGUST | SEPT. j OCT. | NOV. 

« l»» no ItT Ilf Ml MB IBB IBB iff I7B Mi fllMf7 BOB tM MB ttS tit Mt B4B IBI tlOtBB BT4 111 IBB BtB BOB BOB IM Ut Jt4 

ti BO 7 II tl tB t 4 II IB tB B t IB It 10 B IB to CT i 10 17 Bl I B IB ft IB B tB M 10 
* " r ' I "» I * I ■■ I ■ t < 4 I t I ^ I I I I I , I I I 






* JULIA*! OATES 

Source: U.S, Weekly Crop and Weather Bulletin 
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The experimental findings suggest that the driver of discrimination 
is not spectral shape* rather the relative difference of amplitude of the 
spectral mean values* which occurs when the crop species are observed at 
different growth stages. 

This does not connote the absence of all Information in the quantitative 
differences between diverse spectral bands. It simply says that the infor- 
mation contained within the inter-band variations is relatively less impor- 
tant than that Inherent in the differences between the intensities of the 
mean reflectances. Strt^d in other words* if the spectral reflectance 
shapes of various cro^ or of a single crop at different growth stages* 
are normalized to a common average reflectance* the ability to discriminate 
between them is limited. Interband differences probably do arise* espe- 
cially if the entire region of interest is Imaged simultaneously, from 
momentary differences of leaf orientations, different fraction uf under- 
lying ground covered by the leaf canopy* differences in the reflecting 
properties of the underlying ground, and similar phenomena depending on 
geometry or upon electromagnetic wave Impedances and transmissivities pre- 
sented by the pattern of reflecting objects. 

Nevertheless* the experimental data appears to indicate that at- 
tempts at exploiting differentials of spectral shapes in and by themselves 
is akin to attempting to distinguish targets from clutter using amplitude 
patterns alone. The resulting accuracies are limited. 

The presence of the statistical fluctuations* typical of clutter 
phenomena, further reduces the discrimination capability. The phenomeno- 
logy appears formally equivalent to the glint effect observed in radar. 
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The situation suggests* as a working hypothesis* a model equivalent 
to the blstatlc radar mechanism. The sun provides the Illuminator* the 
illuminated crop represents a "scintillating'* target* and the sensor is 
the receiver. In this model* the discrimination error Is analogous to 
a false alarm rate. 

The hypothesis leads to the development of a generalized theory of 
discrimination* presented In Chapter 5 following. 
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5.0 THEORY OF DISCRIMIMATION FROM RADIOMETRIC DATA 


5,1 Simplified Theory 

The objective of agricultural acreage mensuration Is to determine the 
proportion of the crop or crops sought relative to the total area under 
analysis* Independently of where each crop species is located geographi- 
cally. A first approach to the theory can be made by Invoking the fol- 
lowing assumptions: 

• Two species only; normal distribution of reflectances 

• One radiometric band only 

e Area under study sufficiently large so that mensuration errors 
Induced by the finite size of the resolution element (pixel) 
are negligible 

It is shown In Appendix E that* under these assumptions* the error 
attributable to discrimination In estimating the proportion of the crop 
sought Is given by: 

e % - — /p a + (1-p) a (5-1) 

py^ <Pa - 

where : 

p <= proportion of crop A to total area 

e^% <= percent error In p attributable to discrimination 

k = statistical confidence multiplier 

» total number of pixels In scene (species A pixels + 
species B pixels. 


n 
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p^f <■ aeans of the relative reflectances of species A» B 
respectively. 

a. » On standard deviations of the fluctuations of p.« p_ 
A B A B 

respectively. 

Expression (5-1) can also be written: 


100 k h, 




/ 




P <i-p) (ir) 


^ (1 - -^) 


(5-2) 


where: 


h* * — “ coefficient of variation of p, 
A ‘'A 


'B 


h„ •= — •= coefficient of variation of p 

B Pn 


B 


B 


all other symbols are as defined previously 


i^B 

— ' * is the dont/t/l&t Kado between the reflectances of the 
^A 

two crops 


5.2 Significance 


There is currently no definitive experimental indication that the 
coefficients of variation differ systematically and reliably among crops. 
Until more comprehensive data become available, it is assumed that 
h^ = h^ = h. 

Expression (5-2) can thus be modified by posing h^ » h^ ■= h. This 
assumption converts (5-2) into: 
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p ^ <1 - * 


Expllcltatlng n from (5-3) yields: 


(5-3) 


n 



(e/,)^ (1 



Ip + (1-p) (~) 1 


(5-4) 


Let ue survey briefly the implications of this expression. 

p is a characteristic of the target crop, in fact it is the quantity 
to be measured, p does vary from year to year: in most practical cases, 

the variation is relatively small. More significantly, p lies beyond the 
observer's control. He can however render the computation of p more facile 
li he. pfvbjy to ceAtcUn pAlon. tn^o^ation, 

e^%, the allowable one-sigma error, is established by the user re- 
quirements. It is also beyond the observer's control. 


Pa* Pn ^ determined by the characteristics of the crops being 
A B 

inventoried. Control over p., p_ la possible by selecting the timing of 

A B 

the observation. 

n represents the nxmiber of pixels required to be observed. It is 
analogous to the "manber of looks" in radar theory. It is controllable , 
by varying the extent of the surface area being observed. 

n is a driver of the characteristics and cost-effectiveness of the 
Data System, because: 
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i' ^ 


1. The Data System's load is directly proportional to the nuinber of 
pixels to be observed. 

2. The costs of the data relay, preprocessing and Information extrac 
tlon portions of the Data System grow in proportion to the data 
load. 


A feel for the relationships between parameters is available from 
Figure 5-1. It plots the total area required to be observed and measured 

Pfi 

versus the contrast ratio — , for the typical conditions: error desired 

2%; confidence one sigma (k«l}; p«0.15 (a typical realistic value for 
densely cropped areas) ; and a linear pixel dimension of 30 meters (pro- 
posed for LMIDSAT D). As shown in Chapter 4, the values of h used as 
parameters, e.g. 0.1 and 0.3, are typical of minimum and maximum values 
competing to realistic agricultural scenes. 


Figure 5-1 can be converted to the characteristics of LANDSAT A, B, 
C, by increasing the ordinates by the square of the ratio of the respec- 
tive resolutions, i.e. approximately fivefold. 


B 

A feel for typical practical values of the contrast ratio — can be 

^A 

obtained from Figure 5-2, derived from the experimental data summarized in 

Chapter 4. It shows that under conservative system regimes, taking into 

account compression effects Induced by the atmosphere, one should not 

^B 

assume design goals for — significantly lower than 0.8; and that 0.95 

^A 

represents an approximate upper bound, beyond which the area which must 
be measured generally transcends practically attainable limits. 


FIGURE S-1 

AREA WHICH MUST BE OBSERVED AS A FUNCTION OF SCENE PARAMETERS 


DISCRLMINATION EPUIOR, e = 2 \ 

d 

PROPORTION p » 0.15 




FIGURE 5-2 



Barley vs. Wheat Oats vs. Wheat Alfalfa vs. Wheat Pasture vs. Wheat 


Note: All wheat is spring wheat 
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5.3 Selection of the Best Band 

The simplified theory deals with a single band. Given a choice among 
various bands. It comes naturally to ask which band Is Indeed the best and 
how to select It. 

The most efficient system is obviously the one which provides the best 
performance at minimum cost. From the aforegoing, this Is tantamount to 
require the simultaneous minimization of the error and of the number of 
pixels which must be sampled. From expression (5-2), It is clear that 
the product e^% can be taken as representing a figure of merit of the 
Data System. Thus, the "best" single band Is the one which, for specified 
conditions of observation, minimizes the figure of merit e^% /n. This is 
equivalent, from (5-2), to minimizing the expression; 


h / p 2 h 2 

^ /p + (1-p) (^) (^) (5-5) 

^ ^ A 

where the constant K subsumes all parameters not affected by the choice 
of band. Under the previously Introduced assumptions h^ » h^ « h, 

(5-5) becomes! 





r. 


'V + (1-p) (-^) 


2 


(5-6) 
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In other words* Che "best" band consistently exhibits combinations 
of p^, pg, h^, hg which minimize expression (5-5) — or (5-6) if one does 
not have available the separate values of h^« hg. 

Tests were run to determine the best band or bands from the available 
experimental data. To this effect* pairs of crops were chosen whose ground 
cover reaches substantial values* in excess of 50% during the season. The 
reflectance values of the fine-grained spectra were Integrated within the 
wavelength intervals of the first five TM bands. 

From (5-6) * the number of pixels required to be observed at each of 
the available temporal dates were calculated. Typical experimental re- 
sults are shown in Table 5-1 for h 0.3* error e^ » 2%* p « 0.15* crops 
alfalfa versus barley* Hand County data. It is clear Chat the bands 
which require the lowest number of pixels* l.e. which minimize the "cost 
function*" are the best. 


TABLE 5-1 

number of Pixns (cost function) required to be OBSERVEO within toe first fIVE TM BANDS 





DATES* 

(1977) 




OATES. 

(1977) 


TM BAND 

WAVaENGTH. urn 

4/21 

5/10 

6/1 

6/16 


4/21 

5/10 

6/1 

6/16 

1 

0.45 • 0.52 

at 

2.9 X 10® 

2.9 X 10® 

2.3 X 10® 


- 

34 

345 

2.7 

2 

0.52 - 0.60 

8.3 X 10® 

1.7 X 10* 

4.6 X 10^ 

4.9 X 10® 


99 

202 

5»480 

S.8 

3 

0.63 > 0.69 

1.6 X 10® 

1.9 X 10® 

1.0 X 10® 

6 X 10® 


190 

226 

119 

7,1 

4 

0.76 ' 0.90 

3.2 X 10® 

2.6 X 10® 

3.2 X 10® 

B.4 X 10® 


51 

3.1 

38 

Z] 

5 

1.55 - 1.75 

4 4.5 X 10^ 

2.7 X 10® 

4.0 X 10® 

2.5 X 10® , 

t5t360 

32 

46 

30 


NUMBER OF PIXELS COST FUNCTION NORMALIZED TO 

REQUIRED FOR DISCRIMNATIOK LOWEST VALUE 



The example of Table 5-1 shows that TM Band 4 Is the be8t '((o4 

C06C ^eZe.cte.d. 

In general, the best band is case-dependent. Various tests run for 
diverse crops at 50% or better ground cover resulted in the orientative 
rankings of the discriminating power of the TM bands depicted in Table 5-2. 


TABLt 5-2 

RANKING OF TO 8AHDS POR DlSCRimHMlLITY BETW£EW STAPLE CHOPS 


TM BAND 

WAVELENGTH, MrCRONS 

APPROXIMATE 
RELATIVE 
FIGURE OF MERIT 

1 

0.4S • 0.S2 

O.T 

2 

0.52 - 0.50 

0.1 

3 

0.63 - 0.69 

0.15 

4 

0.76 0.90 


5 

1.55 - 1.7S 

0.3 


5.4 Effect of Multiple Bands 

Rigorous extension of the theory to multiple bands was not attempted 
in this effort: because the objective was to determine the drivers of 

Data System performance rather than to develop discrimination theory. 

Some general considerations on the effect of multiple bands can however 
be made at this time: 

1. The addition of other bands to the best band will generally improve 
the discriminant power of the system; provided, that the slgnal-to- 
fluctuation ratio this added to the best band is not so low as to 
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yield a decreased overall signal-to-fluctuatlon ratio for the 
multiple bands. 

2. The addition of bands la beneficial - subject to the caveat ex- 
pressed above - even if their spectral amplitudes are completely 
correlated (cf. Table D-2, Appendix D). This is because adding 
bands provides additional samples. Thus, if all bands possess 
the same discriminating power* the upper bound of the improvement 
is given by the sampling theorem* l.e. it will equal /T* where 
b is the number of bands combined. In practice* as indicated by 
Table 5-2* the benefit will be lower than this upper bound* if 
the best band has been selected Initially. 

In conclusion* all other conditions remaining invariant* multlspectral 
discrimination can be expected to lower the number of pixels and corres- 
ponding total area required to be observed (cf. Figure 5-1)* but not by 
a strikingly large amount. 
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6.0 CONSEQUENCES OF THE THEORY ON THE CHARACTERISTICS OF THE DATA SYSTEM 


6.1 Impact on the Frequency of Observation 


The developments of previous Chapter 5 lead to the following 
conclusions regarding accurate discrimination: 


It is not sufficient to image agricultural scenes "when all 
crops look good." They must be Imaged when the statistical 
"distances" between their probability distributions reach suf- 
ficiently large values. The requires selecting moments of ob- 

servation for which — and h. . h„ are, simultaneously, suf- 

p A :b 

flclently small. 


2. The above conditions do not occur for all crops simultaneously: 
rather, they occur generally pairwise. The net result is that 
quality discrimination must be achieved by analyzing several 
optimally positioned temporal "cuts." 


The typical situation is depicted in Figure 6-1. The time Intervals 
during which the "statistical dlstr.nces" are widest bear the generalized 
designation of "temporal windows", or simply "windows." 


The duration of a window is not a unlvocal quantity. It is a function 

of : 

1, The system's specified figure of merit, e% The more stringent 

the figure of merit i.e., the higher the desired accuracy and 

simultaneously the lower the number of pixels one wishes to pro- 
cess — the smaller must be the value of function (5-5). Conse- 
quently the shorter will be the admissible window. In extreme 


REFLECTANCE FACTOR, PERCENT 


FIGURE 6-1 


CONCEPT OF STATISTICAL WINDOW 


DATA FROM HAND CO., N.D. 
SPECTRAL BAND: 0.76-0.90 


son 


40 ] 


30 ^ 


20 ^ 


10 -^ 


ALFALFA (at 80% GROUND COVER) 



SPRING 
WHEAT 
/ (AT 30% 

GROUND COVER) 


APRIL 


MAY 


JUNE 


JULY 


rn 

o 


ONE-SIGMA INTERVAL OF FLUCTUATIONS 



cases* l.e. If the specification of the figure of merit Is made 
too stringent* the window may vanish altogether. 


The system's sensitivity threshold ^^e previous deriva- 

tions disregard exogenous fluctuations acting within the scene. 
In practice* such fluctuations exist. Thus it is necessary in 
real situations to allow a "margin of safety" on the threshold 
value of In other words* the actual ~ must be chosen some- 

what more Stringently* or lower* than its theoretically accept- 
able value. 


The principal exogenous fluctuations stem from three sources: 

o The granularity of the sensor*' caused by quantization 
errors* calibration imperfections, instrumental errors. 

o The variability of the atmospheric absorption across the 
scene. 

o The variability of the path radiance across the scene. 

As regards the sensor, a reasonable estimate* derived from MSS 
performance, is to assume that the theoretical maxlmtim number 
of levels is reduced by a factor of two. Thus* for example the 
Thematic Mapper* which possesses a theoretical maximum of 256 
levels* probably contains only approximately 100 unlvocal levels 
at a false alarm ratio of order 1%. 

Because experimental results shew that the reflectance fluctua- 
tions of crop targets* also known as "scene noise*" are consi- 
derably larger than this value* sensor fluctuations are negli- 
gible to the effects of performing agricultural inventories by 
radiometric remote sensing from OSTA space platforms. 

Pb 

Atmospheric transmissivity fluctuations affect ( — are those 

Pa 

which occur between diverse points of the scene. Variations 
which affect the entire scene uniformly affect both and 
by the same ratio and thus Impact only slightly the ratio 
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Information on the lntra-*scene fluctuations is limited. Invest!- 

gatlons by EGOsystems* using data from the US Turbidity Network. 

indicate that fluctuations between separate points of the scene 

can reach peak excursions as high as 20%, rms values of as much 

as 10%. Figure 6-2, The contribution of these fluctuations to 

the uncertainty In ( — . depends upon the distribution of the 

p* tn 

species A, B across tne scene. If A and B are adjacent to each 

other everywhere (locally compact), the fluctuations will have 
^ R 

no effect on . If however the locations of A and B are 

P* th 

disaggregated (ror example, as a limiting case, all the A's are 
In one place, all the B's In another, distant from the first), 
then and will suffer different random variations. Pending 
the availability of supplementary information on atmospheric 
fluctuations. It is not unreasonable to assume that p., are 
both affected by uncorrelated random fluctuations of 3%: their 

combined rms fluctuation thus would be of order 4%, one sigma, 

* The third factor, variability of the path radiance. Is even less 

known. It Is not Included In this calculation. 

3. The amount of prior Information available. If the phenologic 
progress for the year under study is poorly known, the average 
phenologic progressions are computed from historical statistics. 
Figure 6-3 Illustrates the procedure for determining the temporal 
fluctuations of the phenologic stages from historical data, In 
this case from [ 10] . it Is noted that the yearly progressions 
of mean reflectances are shown as reaching the same peak values. 
This Is not necessarily true In nature: it is an assumption 

dictated by the scarcity of experimental data. In practice. It 
can be expected that the peak values will vary from year to year, 
in response to differing sltuatfions of ground cover, plant stress 
and similar. The computation of the temporal excursions Is how- 
ever affected but little by fluctuations of reflectance levels. 


11 V-.V'.-'rir' — ... ■ , ... ............ .. '.'•.•j.i. 


Differ 


FIGURE 6-2 



. -Hi. 



FIGURE 6-3 


PROCEDURE FOR CONSTRUCTING PHENOLOGICAL STATISTICS 

RELATIVE 

REFLECTANCE XM BAND 4 



+ SEVEN ADDITIONAL PAST YEARS 


h* = 
ta B 
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HEADING 

TASSELING 

PODDING 

EMERGING 



WHEAT 
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SOYBEANS 


0.6 



10-YEAR AVERAGE 
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4< The presence of cloud cover. The random occurrence of clouds 
during Che window intervals can be calculated or simulated with 
existing cloud models. This is further discussed in Chapter 7. 

Figure 6-4 depicts typical windows. Note the effect upon window 
length of varying the sensitivity threshold of the system. The 
windows discussed previously are for cloudless conditions. Figure 6-5 
depicts the required frequency of observation for selected regions of 
the world when cloud cover Is taken Into account for various levels of 
threshold sensitivity. Chapter 7 Illustrates how the windows portrayed 
in Figure 6-4 for selected regions were expanded worldwide by simulation. 
The results were analogous to those portrayed In Figure 6-5: recurrence 

times of order 5 days are needed to achieve coverages commensurate with 
96% to 98% completeness of crop acreage mensuration. 

Chapter 8 discusses ways and means of achieving such low recurrence 
rates with minimal cost impact upon the system. 


6.2 Impact on the Extent of the Area which must be Observed 

The preceding theory shows that it is not necessary to observe* sample 
and process the entire surface of the region under study. If the times of 
observation are selected correctly, sampling of considerably smaller por- 
tions is sufficient. Let us examine this in more detail. 

6*2. 1 The Effect of Bounding and the Role of Auxiliary Information 

The sensitivity of the area to be measured to its crop proportion is 
depicted in Figure 6-6, which also presents for comparison the proper' ions 


FIGURE C-4 

STATISTICAL (ONE SIGMA) TEMPORAL WINDOW ENVELOPES OF REFLECTANCE DIFFERENCE 

USSR, UKRAINE 



PMaicM. niRESHOiD uetm lj||}{]] 
(SENSM > ATHOSPHCRE) IIIHH 



FIGURE 6-5 

REQUIRED RECURRENCE OF OBSERVATION TO ACHIEVE 99% PROBABILITY OF VIEWING FOR SPECIFIED 
CLOUD COVER AS A FUNCTION OF SENSOR THRESHOLD SENSITIVITY Ap 



aOUD COVER PR 06 AB 1 UTY, % 


AREA TO BE OBSERVED 
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competing to selected areas of the world. It can be seen that the average 
proportions of crops are relatively smallj calling for the sampling of 
rather extensive areas: but that the area to be observed decreases dras- 

tically as the crop proportion Increases. 

This leads to the concept of "bounding,’* l.e. eliminating from the 
observed scene as many as possible areas known not to contain the crop 
or crops sought. The concept of bounding is obvious for such surface fea- 
tures as lakes, deserts, dense urban settlements. For other than these 
obvious features, proper bounding requires the availability of prior In- 
formation, l.e. what areas under observations do or do not contain the 
crops sought. After bounding Is accomplished to the best extent possible 
with the information available, stratification into diverse levels of crop 
I proportions is effective in further reducing the area to be observed. 

Same as with bounding, proper stratification also requires prior Informa- 
tion. 

6.2.2 Sampling 

If the moments of observation are correctly chosen, the area required 
to be observed will be significantly less than the net scene area still 
available, even after bounding and stratification. In other words. It 
will not be necessary to observe all of the residual area: observation 

of a smaller fraction generally suffices to yield the desired accuracies. 
(5-4) allows computation of the area required to be observed: however, 

it does not provide information as to what geographical locales within 
^ the scene should be observed. 
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In the absence of other prior information, the distribution of the 
locales to be observed Is determined by statistical sampling. The stan- 
dard procedure ii; to choose the nwnber of samples, or segments, from the 
well-known formulation: 



\^ere: 

e «■ sampling error desired In the measurement of p 
N -> number of sample segments 

d « fluctuation of proportion within the sampled area 
Q a mean value of the proportion measured from the N samples. 

In practice, o, or, equivalently, the coefficients of variation ~ 
of agricultural regions are not published nor available. What is 
available Is the historical sequence of the region's proportions of tar- 
get crops p. If the successive p's do not fluctuate too widely, one can 
assume their temporal average as representative of the region; other- 
wise, In absence of other information. It Is prudent to assume the 
lowest value of p experienced over the last several years. 

This historical value of p can be utilized in place of the un- 
available ~ by assuming that the population of samples to be chosen 
obeys a Bernoulllan distribution, wherein each sample, or segment, is 
counted as zero If it contains a proportion of the target crop less 
than p: and as 1, If more. With this assumption, the standard de- 
viation 5 is given by: 
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d « p. (1-p) (6-2) 

l.N 

E 

1 

The mean p of the sample is , which, for N sufficiently large, 

tends to p. The mean proportion containing the species sought Is then 
Npt and expression (6-1) becomes: 


^ . , ■^ E. a . 71} . /'S -k-. (6-3) 

s Np ^ S 

Expressing e as a percent error* and Introducing the confidence multi- 

3 

pller k, (6-3) becomes: 

e “ 100 k / (6-4) 
s Np 

This expression is valid if Npk » 1, l.e. to the extent to which 
the Bemoulllan approximation is close to the Gaussian distribution, and 
if N « total number of sample segments extractable from the observed 
scene. In practice, the errors are negligible at the 1% as long as 
Npk ^ 30, and if N/ (Parent population N) _< 2%, 

The formulation above applies to unstratified sampling. The effect 
of stratification is to Increase the "effective" p by parceling the 
samples according to the p - density of the various strata within the 
scene. 

Expression (6-4), or its originating form (6-1), do not in and by 
themselves answer the question of how large to make each Individual 
segment. As the segment's surface is made progressively larger, it 
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eventually will encompass the entire area to Inventoried. The number 
of segments becomes unity s yet the error becomes zero; thus (6-4) 
breaks down in this extreme case. A quantitative formulation requires 
knowledge of the distribution of p within the region. In Its absence, 
(6-4) still holds for conditions wherein the total area spmpled Is *'suf- 
ficlently" small with respect to the total area being Inventoried, and 
the number of segments is "sufficiently" large. 

From (6-4), the number of sample segments can be calculated: 


N « — (6-5) 

(e/)^ P 

We have seen that the total number of pixels required for discrimination 
is given by expression (5-4). The question is, how many pixels to al- 
locate to each segment, l.e. what Is the optimum area of each segment. 

If there are no other error sources, and if the errors Induced by 
sampling and those Induced by discrimination are uncorrelated, then the 
total error ia< given by: 


®t ” (6-6) 

where e is the Inventory error specified by the user, e Is the sampling 
V S 

error, and e^ Is the discrimination error. 

Expression (6-6) Indicates that there is a choice in the selection 
of tl, number of sample segments, and n^, total number of pixels to be 


1 
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observed; as long as neither of the corresponding errors e^, e^ equals 
the total error e^t and subject further to the obvious constraint N < n. 

For exaiiq>let suppose that one vishes to minimise n^, which is pro- 
portional to the data system costs. To see what happens* let us re-write 
(6-6) using expressions (6-4) and (5-3): 



2 

P 


a - 


tP + a-P) (^)'] + ^ 


Using for convenience: 


2 ,, ^Bv2 

p (1--) 


tp+(l-p)(-^)] =A 


P 

and substituting in (6-7) : 



whence : 



(6-7) 


(6-8) 


To find the minimum n^> since e^* A* B are constants, set the condition: 


+1' 
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AB 


N 


(e ^ 

' t r 


(6-9) 


Condition (6-9) is satisfied when the denominator tends to infinity* 
i.e. vdien N -*■ <». 


Since* by definition* N cannot exceed n^* the minimum value of n^ 
is achieved when n^ ■ N. From (6-8)* e^ becomes: 



e 


t 



✓ A + B 


In this case* each pixel acts both as a discrimination element and as a 
sample segment. The total number of required pixels will be given by: 



n 


A + B ,„4 

ST- « 10 


A 4- B 
(e^%)^ 


(6-10) 


and will be somehwat larger than depicted in Figure 5-1. In practice, 
this mode of operation may be difficult to achieve by orbital sensors 
of the not too distant future* for two reasons: 1) it would entail a 

very rapid scan of the sensor during overpass; 2) it would entail high 
precision of pointing* to Insure correct location of the "pixel segments" 
on the earth's surface*. 


On the other extreme* if one wishes tb minimize the number of sample 
segments N» so as to minimize the sensor's angular excursions* it can be 
observed that N's lower limit is given by (6-5). However* choice of 
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thls value of N entails n ■ because it is tantamount to specifying 
e^ •> 0 in (6-6) and hence in (5-4) . In this extreme case* the area of 
each segment would become infinite. Since this is not possible, the 
upper limit for this case is the already-discussed case of a single 
segment encompassing the entire area to be inventoried. As N is in- 
creased from this lower limiting value, n decreases monotonlcally until 
it reaches the value given by (6-10) . 

Optimization of the relative weights to be placed on e^ and e^ 
requires the formulation of an additional relationship between N and n. 
Such a relationship could be written, for example, if onw knew the im- 
provement of sampling accuracy achievable by increasing the area sub- 
tended by each segment. This type of relationship would be a function 
of the characteristics of the observed region. 

In the absence of such knowledge, one can assume that e and e. 

s d 

possess equal weights in the total error budget: 



With this hypothesis, the total number of pixels required, n^, is, 
from (5-3) : 


4 2 2 

2 X 10^ k"^ 

2 2 2 
(1 - 


(p + (1 - p) (^) ] 


( 6 - 11 ) 


and the number of samples, N, is from (6-5): 
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N. 2 X loSVcV- P) 


( 6 - 12 ) 


p 


Consequently* the number of pixels per sample segment* n * becomes: 

8 


n 


2 


n 


8 N 






P <1 - 3 


^1 - p 


+ (~) 1 


(6-13) 


Figure 6-7 shows n as a function of p for values of the parameters 
Pb ° 

e^t h* corresponding to practical situations* The strong effect of 
increasing p on reducing the sample areas ia clearly evident. The 
corresponding mmiber of segments* N> is shown as a function of p in 
Figure 6-8, 


As already mentioned* the computations above hold for unstratlfled 
sampling. Stratification generally reduces the number of samples re- 
quired. USDA* for example* estimates (personal communications from 
Mr. William Wlgton of USDA/eBCS) that the number of samples required 
to inventory U.S. wheat at a 2% sampling error is of order 1*000 with 
optimum stratification. This number would Increase If one wished simul- 
taneously to Inventory additional crops. 


The allocation of stratified samples is a function of the homo- 
geneity of the areal distribution of acreage. To assess the latter, 
a map of the region's agricultural acreage Is required. Because this 
varies with time* sampling schemes periodically require the availability 
of full frames, at intervals commensurate with the dynamics of the ob- 
served region: generally of order several years. 





TOTAL AREA SANPUO, llicurit 
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6.2.3 Mensuration 

If the total area to be Inventoried Is relatively small* an addi- 
tional source of error* well known In surveying practice* arises due 
to the finite precision of the ireasurement devices employed. In the 
case at hand, the precision Is of order one pixel. The corresponding 
mensuration error Is: 


e * s IXi-dL . ( 6 _ 14 ) 

“ ^ /n p «■ 

Where: 

e^% “ mensuration error, percent cT crop area Inventoried 
Y form factor of sample (« 1 for square samples) 

3 = mensuration figure of merit (t 1 for photointerpreter) 
r = linear dimension of pixel, meters 
k “ statistical confidence multiplier (= 1 for lo) 

A = area of each segment, hectares 
p = proportion of crop area to total area 
N = number of sample segments 

n = total number of pixels in all sample segments 

6 = area of each pixel, hectares (0,09 for 30 meter 
LANDSAT D pixels, :0.44 for LANDSAT C pixels). 


The mensuration error, depicted In Figure 6-9 for practically en- 
countered parameters, must be taken Into account when it becomes a sub- 
stantial fraction of the allowable total error e^. This can happen if 
the total area to be sampled is small — brought about for example, if 
p is very high and/or the contrast ratio — is very low. 
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In this case the total error budget must be allocated between the 
three sources of error* i.e. sampling* discrimination* mensuration; 

2 2^2^2 1CX 

e. •• e + e . + e (6-15) 

t s d m 

In the absence of other prior knowledge* the error budgets are allocated 
equally. The simple manipulations result in expressions similar to 
(6-11), (6-12), (6-13), with the additional term due to mensuration 
error included. Figure 6-10 depicts the three sources of statistical 
error arising in proportion estimation, and the corresponding expres- 
sions for each and for the total error. 

6.2.4 Advantages and Practical Constraints of Sampling 

Sampling can significantly reduce the systems data load: specifi- 

cally* in spacecraft-to-ground data relay, in preprocessing, and in the 
information extraction function. 

The reduction theoretically achievable by sampling vis-a-vis 
wall-to-wall imaging is shown in Figure 6-11. In the near future, 
technological constraints limit its full realization. For example, 
the precision of the platform's stabilization system impacts the ac- 
curacy with which the sample segment can be located on the earth's 
surface. The typical stabilization accuracy of the Multi Mission Mo- 
dular Spacecraft (MMS), planned standard support for LANDGAT-D, is 
±0.25 arc min rms. From an orbital altitude of order 700 Km, the cor- 
responding ground smear is of order ±140 meters, 99% of the times. 
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REDUCTION OF DATA LOAD THEORETICALLY ACHIEVABLE BY SAMPLING 
VERSUS WALL-TO-WALL COVERAGE 
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Sampling cannot be used to the total exclusion of all wall>to-wall 
data« The latter Is still required periodically to update the strati- 
fication method. 

6.3 Role of Auxiliary Information 

The preceding developments illustrate how crop discrlrjiinatlon improves 
if prior available information is exploited in conjunction with radiometric 
data. In particular* the availability of prior information appears essen- 
tial to the crop function. 

Limited experimental results indicate the possible appearance of 
narrow band spectral features peculiar to particular species at some 
phenologic stages. The reliable presence of such "indicators" would 
allow identification without prior information. So far, the only such 
effect known to the authors to have been observed is the "red shift" 
of wheat, at the pre-flowering stage, found by S. Ungar of GISS [15]. 

The value of prior information to the discrimination process is 
greater, the less redundant its contents with the radiometric informa- 
tion supplied by the sensor. -The much used procedure of deriving basic 
spectral data from training sites is thus valuable only if natural crop 
spectra differ significantly one from the other. According to the model 
•previously described, this occurs primarily during selected time periods 
corresponding to the phenologic windows. It is in these periods that 
training site information has maximum value. 
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Errors of mensuration Induced by the coarseness of the pixel dimen- 
sion can be reduced by appropriately combining the low-resolution radio- 
metric data with high-resolution aerophotographlc data. The latter 
serve to delineate geometrical boundaries » the former to Identify their 
internal contents. 


Table 6-1 presents the types of auxiliary information which most 
significantly contribute to the accuracy of crop inventories. 


TABLE 6-1 

pRiHcr ;.; auxiliary data peculiar to heciom of ihtctest 

CROP ACREAGE INVENTORV 

BASaiNE DATA 

0 STATISTICAL DATA1 PRODUCTION. ACREAGE. PROPORTIONS OF CROPS; FIELD 
SUE STATISnCS. 

0 PHENaOGIC CALENDAR ITS STATISTICAL FLUCTUATION 
o REPRESENTATIVE OR INDICATIVE SPECTRA 

o HETEOSTATISTICS, RAINF/U.L. INSOLATION, TEMPERATURE, EVAPOTRANSPIRATION 

o aOUD COVER STATISTICS 

o TOPO HAPS, VARIOUS SCALES 

o UTEST CONVENTIONAL AEROPHOTOS 

o LATEST LAND USE HAPS 

<1 AGRICULTURAL PRACTICES, AGRt ECONOMIC DATA 

0 STATISTICS ON PRINCIPAL CAUSES OF CROP LOSSES 

o DESCRIPTION AND DATA ON CONVENTIONAL (GROUND-BASED) INVENTORY METHODS 
0 CADASTRAL HAPS 
DERIVED DATA 

o AGROMET INDICES 
o ACREAGE. Viao, PRODUCTION TRENDS 
0 SIGNATURE EXTENSION INDICES 
SOFTWARE h 

o CONVENTIONAL PROGRAMS NODULES AND PACKAGES 
o SPECIAL PROGRVI CLEARINGHOUSE 
DATA-REIATED SERVICES 

o RECOtflENDATIOH OF OPTIHUH SCENE TIMING 
o RECOmENDATION OF TUNING PARAMETERS FOR SPACE-DERIVED DATA 


7.0 SIMULATION SCHEME AND RESULTS 


The approach orglnally conceived was to test the theory thus far 
presented on a worldwide basis, including the effects of sampling, 
phenologlc cycles and cloud cover, by simulation on the Marshall 
Space Flight Center's Data System Dynamic Simulator (DSDS). The 
DSDS concept embodies these capabilities. Including capability for 
radiometric readout from imagery, magnetic tapes, or simulated scenes. 

At the time of this simulation, the DSDS had not reached its Phase 
III development stage. Hence the simulation was tailored to accept, 
instead of the full complement of inputs in their orglnal format, 
the following: 1) Synthetic sampling and discrimination derived from 

theory; 2) Synthetic model of the reflectance excursions induced by 
phenology, derived from empirical data; 3) Pre'-determlned geographic 
locations of the principal crops, derived from historical statistics; 

4) Synthetic cloud cover statistics, derived from existing models. 

The simulation scheme was elaborated jointly with NASA Marshall and 
with the General Electric Company's Huntsville Operations, who ran the 
simulation. General Electric's results are contained in separate report 
[8}. This chapter synopslzes bow the simulation inputs were structured, 
and reports key results. The procedure adopted was the following: 

1. The principal comtrles and. crops were determined from the 
statistical material presented in Chapter 2. Twenty-two 
countries were selected. Pre-simulation analyses showed 
that adequate results could be obtained by simulating only 
the crops Winter Wheat, Spring Wheat, Corn, Soybeans, and 
Rice . 
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2. For each selected country, phonologic calendars were obtained 
from I9J. 

3. Because not all crop calendars provided complete Information 
on all phenologlc stages, they were adjusted against the 
data of the SRS Crop Weather Bulletin [10]. From these, 

U.S. phenologlc statistics were compiled following the 
procedure described in Chapter 6, Illustrated In Figure 
6-3. The mean and standard deviations of each phenologlc 
stage's time of occurrence was computed. 

Most foreign data were incomplete, e.g. provided only the 
mean dates of occurrence of a few phenologlc stages. The 
phenostages were assumed to possess the same standard 
deviations as confuted for the U.S, Subsequent work by 
Ecosystems [11] Indicates the validity of this procedure 
for countries located In the temperature zones, and for 
tropical countries subject to regularly recurring weather 
phenomena, e.g. monsoons in India. 

4. The phenologlc phasings differed considerably among the 22 
countries selected. Equivalent phenologlc zones were 
defined as those wherein the phasing between phenostages 
did not exceed 20 days. Five such phenologlc zones were 
sufficient to account for the selected 22 countries 

I 

5. The major countries, e.g. the U.S . , USSR, Brazil, turned 
out to contain more than one phenologlc zone. They, were 
thus subdivided Into sub-zones. Thirty-six regions were 
thus defined, as shown in Table 7-1 [83. 

6. The major areas of cultivation of the five selected crops 

were identified and bounded in each of the 36 regions, from 
[14]. Analysis disclosed almost everywhere the presence of 
two major confusers: 1) permanent pastures, l.e. areas 

where grasslike crops are grown for grazing; and 2) hay-type 
crops, such as alfalfa, which are grass like crops harvested 
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TABLE 7-1 

COUNTRIES, REGIONS, CROPS AND NUMBER OF SAMPLES USED FOR STUDY 


COUNTRY & REGION 



Estimated Nb. of bPER- 

CROPS 

SIMULATION SAMPLES 

ATIONAL SAMPLE SEGMENTS 

Argentina 

W, C 

45 

1500 

Australia 

W 

30 

1000 

Bangl adesh 

R 

30 

1000 

Brazil North 

C 

30 

1000 

Brazil South 

C. S, R 

53 

1767 

Canada 

W 

30 

1000 

China North 

W, C, S, R 

60 


China Central 

W, C, S. R 

60 

2000 

China South 

W, C, S, R 

60 

2000 

Egypt 

C 

30 

1000 

France 

W, C 

45 

1500 

India Punjab 

W, C 

45 

1500 

India Ganges 

W, C, R 

53 

1767 

India Central 

W, C, R 

S3 

1767 

India Bilaspur 

W, C, R 

53 

1767 

India Coastal 

R 

30 

1000 

Indonesia 

R 

30 

1000 

Italy 

W. C. R 

53 

1767 

Japan 

R 

30 

1000 

Mexico 

C 

30 

1000 

Pakistan 

W 

30 

1000 

Romania 

W, C 

45 

1500 

South Africa 

C 

30 

1000 

Philippines 

R 

30 

1000 

Thailand 

R 

30 

1000 

Turkey 

W 

30 

1000 

USA - Region A 

W. C. S 

53 

1500 

USA - Region B 

W, C 

45 

2000 

USA - Region C 

W,C, S. R 

60 

2000 

USA - Region D 

W, C, S,R 

60 

2000 

USSR Latvia 

W, C 

45 

1500 

USSR Ukraine 

W. C 

45 

1500 

USSR Transvolga 

W, C 

45 

1500 

USSR Volga-Ural 

W. C 

45 

1500 

USSR Siberia 

W, C 

45 

1500 

Yugoslavia 

W, C 

45 

1500 

TOTAL 


1553 

51355 


W “ Wheat (Winter & Spring); C = Corn; S = Soybeans; R = Rice 
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periodically for bay, sometimes as much as 6 to 7 times 
yearly. The phenologlc cycle of these two types of 
confusing crops affects the discrimination of the 
primary crops differently. 

Permanent pastures generally reach the stage of maturity, 
l.e. of high reflectance, prior to the emergence of the 
wanted primary crops. It is thus necessary to "mask out" 
the former from the remainder of the scene by employing 
LANDSAT data gaterhed sufficiently early In the season. 

The hayllke crops exhibit early season characteristics 
similar to those of the pastures, yet contribute additional 
information: their frequent cuttings display a characterls- ' 
tic temporal alternation of high and low reflectances. 

For reasons of simulation expediency, the hay] ike crops 
were assimilated within the simulation model to the 
permanent pastures. A corresponding pattern of reflectance 
progression was constructed for each of the 36 regions. 

7. The phenolpglc "windows" were constructed for each of the 36 
regions. The mean reflectances corresponding to each pheno- 
stage were assuned for all regions to equal those observed 
In the available U.S. experimental data. The fluctuations 
around the means were assumed to equal +^30% one sigma, and 
to be uncorrelated among different crops. The unperturbed 
threshold of separation between reflectances was set at 
in 20 years. This means that deviations from the reflectance 
of two crops must .have a joint probability of occurrence of. 
0.05, l.e. individual probabilities of occurrence of (0.05)]g = 
0.22. On the assumption of normal distribution of the 
fluctuations, this yields a threshold fluctuation level of 
0.76 for each of the two crops to.be discriminated or 23% 
of their mean reflectances. The statistical window durations 
were set to the 95% confidence level: l.e. windows shorter 

than chose which were modeled occur on the average less 
frequently than every 20 years. 
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8. An average sensitivity threshold was computed for all win- 
dows, by RSS combinations of : J) 1% for sensor fluctua- 

tions; 2) + 6% for atmospheric fluctuations as described 
in Chpater 6. In order not to overly complicate the 
simulation, the sensxtivity threshold was assumed fixed 
rather than proportional to the reflectance's amplitude: 
its value was selected as 0,02^ or 2% of the relative 
maximum reflectance of 1. T 3 rpical examples of windows 
are depicted in Table 7-2, 


T«LE 7-2 

EXAMPLES OF UINDCU FORMAT USED IN SIKULATION 


COITNTRY 

REGION 

CROPS 

WtNOOU STARTS 
JUUAI4 QATE 

W:>4D0W DURATIONS 
DAYS 

U.S. 

A (IL. 10. OK. 
KS) 

WW-P 

111 

27 


B (NT. SO. NO, 

CO) 

SW-P 

MB 

40 


C (VA, KY. NC, 
TN. SC. CA, 
AL. HI . LO) 

c-s 

151 

17 

USSR 

UKRAINE 

w-s 

115 

46 


WW “ Winter Wheat, SW = Spring Wheat, P » Pasture, C ■ Corn 
S ~ Soybeans 

9. A sampling scheme was constructed within each region. Based 
upon the results derived in Chapter 6, good stratification 
and bounding was assuned: each candidate region was assigned 

1,000 sample segments if it contained a single crop: and 

1,500; 1767; 2000 segments if it contained two three, four crops 
respectively. For expediency of simulation, sub-multiples ot 
these numbers of sample segments were chosen as Inputs to 
DSDS: specifically, 30, 45, 53 and 60 for single, double. 
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trlple and quadruple r:rops respectively. This reduced set 
of segments was selected at random from the larger set,. The 
list of countries* regions* crops and number of sample seg- 
ments contained within each Is shown In Table 7-1. Not there- 
in shoim but Implicit In the previous construction of the 
windows Is the presence of the selected confuser* pasture* 
discussed under point 7* preceding. 

10. Simulation of the complete logic of the cumulation of Informa- 
tion from multiple looks exceeded the constraints placed on 
this effort. In Its place* target acquisition was specified 
as occurring when at least one "look" per phenologic window 
was possible. Whenever such look did not occur* either due 

to excessive shortness of the windows or Impeded by the pre- 
sence of cloud cover* the particular region and crop affected 
was considered as not having been surveyed. 

The "one look" assumption* necessitated by simulation expe- 
diency* is reasonably adequate for the purposes of this pro- 
ject, because: 

a. The specification of a single look per window Is some- 
what lax If the single look is all that is available. 
Ecosystems experience with real scenes shows that 

'generally more frequent looks are needed: an appro- 

ximate average Is of order 1.5 per window. 

b. The requirement of one look per window Is somewhat 
stringent if a sufficient number of earlier or later 
looks are provided wherefrom to reconstruct the state 
of affairs within the missed window. 

11. Cloud cover statistics were Incorporated in the simulation 

by the General Electric Company, in accordance with the available 
worldwide Allied cloud model. The model's statistics are sup- 
plied In terms of percent cloud cover. 
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A more cogent statlstlct l.e.» the "probability of effecting 
target acquisition," requires knowledge of the statistics of 
unobstructed vision within each scene, l.e. of the patterns 
of "holes in the < louds." The bSDS, in its full phase III 
implementation, would have been capable o', deriving these 
statistics from automated analysis of historical LANDSAT 
and/or GOES and NIMBUS scenes. In its absence, it was instead 
decided to incorporate into the simulation two alternate cloud 
cover thresholds: S0% and 90%. When cloud cover exceeded the 

threshold, the acquisition of a crop target was considered 
missed. 


The principal results derived from the simulation for the LANDSAT-D 


configuration are indicated in Table 7-3. 


TABLE 7*3 


PERCENT OF SIMULATE SAMPLE SEGMENTS ACQUIRED BY A SINaE 
ORBITING PUTFORH [6] 

PLATFORM CHARACTERISTICS: 

ORBITAL ALTITUDE: 704 km 

ORBIT INCLINATION: SUN-SYNCHRONOUS 

SNATH WIDTH; IBS km 


COUNTRY/REGION 

PERCENT 


P CLOlfl) COVER b0% 

8 CLOUO COVER m 

ARGENTINA 

93.3 

97,8 

AUSTRALIA 

95.0 

100.0 

BANGLADESH 

84,4 

95.6 

BRAZIL - NORTH 

100.0 

96.7 

BRAZIL - SOUTH 

99.1 

98.1 

CANADA 

96.7 

96,7 

JHINA - NORTH 

B5.8 

88,3 

CHINA - CBITRAL 

76.7 

96,8 

CHINA - SOUTH 

87.5 

90,0 

EGYPT 

100.0 

100.0 

FRANCE 

86.7 

97.8 

INDIA - PUNJAB 

87.8 

98. 9 

INDIA - GANGES 

75.0 

85.4 

INDIA - CENTRAL 

72,3 

90.6 

INDIA - BILASPUR 

82.4 

89,9 

INDIA - COASTAL 

85. 6 

87,8 

INDONESIA 

85.0 

91.1 

ITALY 

86.9 

93.5 

JAPAN 

95.0 

93.3 

MEXICO 

97.2 

98.9 

PAKISTAN 

93.3 

80,0 

ROMANIA 

91.1 

91.1 

S. AFRICA 

90,0 

93.3 

PHILIPPINES 

91,7 

96.7 

THAILAND 

91,7 

96.7 

TURKEY 

93.3 

80,0 

US - A 

60. 0 

60.0 

US - B 

93.0 

93.0 

US - c 

84.6 

87.9 

US - 0 

77,0 

84,7 

USSR - LATVIA 

96.7 

100,0 

USSR - UKRAINE 

96.7 

96,7 

USSR - TRANS-VOLGA 

97.8 » 

95.6 

USSR • VOLGA-URAL 

98.3 

96.7 

USSR - SIBERIA 

94.4 

97.8 
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The Table shows results similar to those presented In previous 
Figure 6**5: namely* a single orbiting platform with the characteris- 

tics of LANDSAT-D -achieves order of 90% of the desired worldwide ob- 
servations. Whereas this result might at first blush appear incom- 
mensurate with the specified accuracy of 96-98% worldwide* It does not 
tell the entire story: for not all countries and regions missed con- 

tribute equally to the world's production of the identified four 
staple crops. A more sophisticated simulation scheme is required to 
determine the weighted Impact of missed observations. 

In any event* modifications to future earth observation systems* 
other than the brute force expedient of multiple spacecraft* appear 
desirable and practically feasible* discussed in the following 
chapter . 
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8.0 IMPACT UPON OSTA*S AGRICULTURAL DATA SYSTEMS OF THE NEXT DECADE 


8.1 State of the Data System 

The reason for being of OSTA Space Systems is to provide information. 
Thus the Data System represents their paramount element » lAich all other 
Subsystems are designed to support. 

To compete with other acreage-measurement systems - e.g. aircraft* 
ground surveys - OSTA Agricultural Data Systems must provide better price/ 
performance. Industry tests the price/performance and market acceptance 
of a given product* prior to its entry to market* by a process known as 
product assurance. A fallout from this study is a preliminary product 
assurance review of OSTA agricultural Information systems. 

The principal design recommendations which affect the overall system’s 
price-performance In Agricultural Applications are recapitulated following. 

8.? System Performance - Accuracy 

To achieve market acceptance* the Data System's accuracy need be' In- 
creased from the currently achieved average level of approximately 80-85%* 
one sigma, to 96-98%, l.e. by a factor of at least three, preferably five. 
Except for the Interference from cloud cover* there are no theoretical 
reasons barring this achievement. Improvements of on-board radiometric 
calibration, better registration of images* and similar sensor-oriented 
procedures* will assist but will not correct the problem. principal 


bottleneck lies In the Information-extraction subsytem: specifically it 
is attributable to the manner in which radiometric information has been 
heretofore exploited. 

Prlnlcpal corrective actions are: 

• Exploitation of the phenologlc windows phenomenon 

• Integration of auxiliary information with the radiometric 
sensing data. 

The first requires significant Increases of the observations 's 
recurrence frequency. The second requires making available suitable 
external information, approximately formatted. 

8.2.1 System Performance-Recurrence Frequency 

The data recurrence frequency must be sufficiently high to cope with 
the finite time duration of the temporal windows, further reduced by cloud 
cover . 

The polar-orbit, sun-synchronous, 185 km swath configuration charac- 
teristic of LANDSAT Liupplies adequate worldwide crop visibility if multiple 
(order of three), suitably phased orbiting platforms are provided. 

Promising avenues to obviate this costly expedient are: 

• Use of microwave sensing in conjunction with optical and near-IR 
sensing 

The technical feasibility of active microwave sensinf from 
orbit was tested in the SEASAT-A spacecraft. 

Limited experimental data in crop reflectances in the microwave 
spectrum indicate strong analogies with their behavior at 
optical wavelengths. Theory supports these findings. 


Increase In the sensor's swath width 


Current thinking holds that nadir viewing of agricultural crops 
is preferable to observations from slant angles. No adequate 
experimental verification of this assumption was found in the 
literature: In fact, limited experimental data [ } point to 

advantages from slant observations. 

The effect of increased swath widths is to lower the repetition 
Interval required of the observations. Its drawback, i.e. re- 
duced cross-track resolution or the need of more sophisticated 
sensors to maintain resolution, can be obviated by replacing 
wall-to-wall with sampled observations. 

Modification of orbital parameters. 

Reductions of recurrence Interval are achievable by modifying the 
platforms orbit. Those are however of lesser magnitude than what 
can be accomplished by the methods previously discussed. 

The largest effect occurs if the orbital inclination is reduced 
so as to cover only the significant agricultural and/or forested 
regions of the world. Figure 8-1 suggests an orbital inclination 
of order 60*’. 

This system modification would entail a departure from sun-syn- 
chronous operation. The corresponding effects upon discrimination 
remain to be tested experimentally. It Is noted that no adequate 
justification for sun-synchronism was found in the literature: nu- 
merous discussions with members of the remote sensing community' in 
dlcate that the specification of this requirement is at present 

founded essentially upon expert opinion rather than upon observed 
fact. 
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8.2.2 System Perfonnance - Role of Auxiliary Information 

It is necessary to combine radiometric imagery with auxiliary 
Information to reach levels of performance comptalble with the user's 
requirements. Most US auxiliary data are avialble from data bases 
managed by Federal Agencies; the residual is obtainable at State and 
Local levels . 

Ihe primary sources of foreign data are agencies at the 
Ministerial level. Albeit with exceptions, foreign data are signifi- 
cantly scantier and of more difficult access than in the US . The 
situation partlculary critical with regard to aerophotography; even 
where available, this most important auxiliary information is so time- 
consuming and laborious to retrieve as to veiy seriously affect the final 
product's timelines and price. 

This situation can be obviated by providing easily accessible 
global photographic base coverage, at appropriate geometric resolutions. 
The, refresh cycle should be commensurate with the period of obsolence 
of the observables: in the case of agriculture, it need be no shorter 
than two or three years. 

A means of implementing this service would be to reduce Che RBV'S 
field of view, coupled with progressively stepped side look. Table 8-1 
shows the times-resolution tradeoffs for complete coverage of th% 
earth surface. 
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TWLt e»i 

Tin atOUIRED FOR COMPLETE BBIf COOEmCE OF THE EABTH‘S 
StlRfACE VERSUS RESOLIITIOW 


RESOLUTIOW. METERS 
40 
20 
10 


C0VERA6E T1HE REOUIBED. YEARS 
0.2 
0.7 
2.85 


It can be seen that resolutions adequate for agricultural appli- 
cations are achievable within practical time constraints. 

i 8.3 System Costs 

The wall-to-wall method of data acquisition is costly on two accounts: 
o It requires very high data relay rates (84 Mblt/sec for LANDSAT D) 
o It requires pre-processing subsystems of high throughput 

The peak throughput requirements are particularly high In certain 
periods of time. During the months from March through August In the 
northern hemisphere, and their counterparts In the southern, the simu- 
lation shows that as much as 500 observations are required daily. 

Preliminary computations £ J Indicate that as much as 50% 
of the LANDSAT D system's running costs Is attributable to the com- 
bination of TDRS channel rental lus cost of pre-processing. 
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As shown in Figure 6-11, reduction of the peak data load by as much 
as 30:1 Is achievable by «&ploylng sampling as the primary data-gathering 
mode. Wall-to-wall data are required to provide !. ase information, but 
only occasionally. 

The presence of reflectance fluctuations Indicates another possible 
avenue of prlce/performance enhancement for agricultural observation 
systems, through relaxation of the sensor's NEA specifications. The 
sensor's sensitivity should be matched to the overall system's error 
budget, of which the "scene fluctuation" is the major contributor. 

Sensor sensitivities of order 6A levels appear adequate. The 
resulting 4.1 reduction of sensor sensitivity can be traded off 
against lower-cost sensor implementation, or It can be exiiloited to 
increase the sensor's geometric resolution. It can be shown theoretically 
that the complete tradeoff would Improve LANDSAT'D resolution to 
approximately 15 meters. 
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APPENDIX A 


THE UM/FAO SCENARIO FOR GLOBAL CROP FORECASTING 


The UN/FAO vlewpolnC can be synthesized by the following logical 
sequence of propositions: 

1. The World's human food supply Is dominated by cereals, which 

account for 53% of the total, followed by 12% tubers, 15% 

hlgh~proteln foods (meat, fish, milk), 20% vegetables, fats, 
sugar combined [ 2 ], 

2 . The global growth of cereal production outstrips population 
growth. Figure A-1 appears to Indicate that no global nutrition 
problem exists [12]. 

3. However, more detailed analysis of the cereal production and popu- 
lation trends of each of the three FAO zones, as in Figure A-2, re- 

veals that the problem Is one of distribution: In the most populous 
Zone C, production Increases barely pace with population growth 

Because Zone C food Intake Is now only approximately 80% of the FAO 
norm. Zone C appears to be hopelessly undernourished. 

4. A deeper step In the analysis, portrayed In Figure A-3, 
indicates that even in Zone C, cereal availability per 

- capita significantly outstrips the consumption, even though the 
apparent per capita consumption Is quite high. 

5. What causes the apparent Inconsistency? Figure A-4 provides the 
explanation. Almost three billion large farm animals, not counting 
upwards of 10 billion heads of poultry, share the earth's surface 
with man. Up until the late fifties, these animals subsisted pri- 
marily on pasture and grass products; in the last decade and a 
half, the composition of feed has been significantly enriched by 
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the so-called ''concentrates," l.e. cereal products: com, soy, 

wheat . 

6. The reason? The growing demand for protein: its efficient and 

cost/effective production requires concentrated fln<Tnal fodder. 

This "technological revolution" has vastly increased the importance 
of cereal products, much beyond the proportions depicted in 
Figure A-1 for human consumption alone. Since meat, milk and ani- 
mal fat production Is now driven by cereals availability. It Is 
safe to estimate that circa 75% of global human food consumption 

is based on cereals. 

7. The cure? Figure A-S indicates that on the average. Zone A*s 
surplus just matches Zone C's deficit. The simple device of 
straightforward transfer would thus, at least on the surface, 
resolve the global food problem. That Is, except for the eco- 
nomic implications, l.e. the matter of who would foot the bill. 

8. The cost of such transfer was computed by the FAO In their re- 
cently compiled World Indicative Plan. As plotted In Figure A-5 
the straightforward transfer of food Including acquisition and 
shipment, would cost $63 billion per year, growing to approxi- 
mately $80 billion In 1985. The cost of financing adequate In- 
frastructures fertilizer plants, farm machinery, education, 
etc. — to enable Zone C to become adequately self-sufficient 
would cost $25 billion yearly growing to approximately $39 bil- 
lion a decade hence. In practice, a complete solution of the 
world's food problem would consist of a mixture of these two 
methods. . In the limit. It would entail a yearly expenditure of 
almost $90 billion yearly currently, growing to exceed $100 bil- 
lion yearly in the mld-1980's. 

9. While these expenditures are- not Impossibly high — the U.S. 
share, for example, of the "maximum plan" would not exceed $36 
billion yearly now, growing to $45 billion In 1985, or approxi- 
mately one-third of our current National Defense outaly — they 
are not likely to be Implemented in the near future. 
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Note that the contributions for agricultural aid to Zone C In all 
forms, from n^l of Zone A and B combined, cover not more than approxi- 
mately 10% of the requirements stated by the FAO. 

The analysis above leads to the conclusion that the primary require- 
ment to alleviate the global food problem Is to manage as best as possible 
the scarce amount of agricultural aid currently flowing from Areas A and B 
Into Area C. 

The requirements ’’universe" as seen from the UN/FAO outlook can be 
divided logically as follows; 

1. Countries in which internal production consistently equals the 
demand, and which are thus self-sufficient nutritionally. 

2. Countries whose Internal production is consistently less than 
the demand, but which are able to import the deficiency through 
other sources of revenue. Typically, the OPEC countries. 

3. Countries with consistent or periodic production-demand deficits 
and which are unable to fill the deficiency by imports, 

4. Countries which are consistent overproducers and therefore 
exporters . 

Countries which exhibit characteristic 1 are not Important to the 
UN/FAO goal, because they already meet it: nor can they directly assist 

in transferring surplus to other diet-deficient Countries. 

Countries exhibiting characteristic 2, l.e. which are consistently 
In crop deficit but which can Import the deficit through other means of 
revenue, are only Important to the extent that their fluctuating import ' 
requirements may Increase or reduce the surplus available to the chroni- 
cally deficient Countries. 






Survey of Countries in Category 3> chronically deficient and unable to 
import the deficit, are of interest to assess the quantity of product they 
are likely to have available in any one year, and thus the extent of the 
current deficiency. This assessment serves as the basis for UN/FAO is 
planning of the amount of food aid to be provided. Aside from political 
considerations, the Interest is proportional to the extent of the deficit 
and the number of persons affected. 

Surveys of Countries in Category 4, consistent overproducers, are of 
primary interest in proportion to the extent of discretionary surplus they 
historically produce. The assessment of their production serves as the 
basis for planning the most appropriate distribution of food aid. 

It is worth noting that the UN/FAO scenario presented above is 
idealized: in current practice, the food aid funds at the disposal of 

the UN are insufficient to cope with the demand. Also, AlD-dlspenslng 
Countries do hot act in harmony and nor for the solely humanitarian pur- 
pose of providing their surplus where most needed': rather, political 

considerations tend to channel food aid where each donor deems it to be 
most opportune. 

In spite of these departures from ideal perfection, the UN/FAO 
scenario is of value as an approximate model to guide the equalization 
between deficits and surpluses. 

Table A-1 dep:'.cts the aggregate distribution of human consumption in 
diet-deficient regions. It shows the presence of the staples already iden- 
tified in the analysis of the U.S. balance of trade situation, plus the 
emergence of rice as a new major staple. 


A- 10 




TflTf r ^'- 


/' 


1 


Table A*>1 also Indicates the 6-year cumulative import-export trends 
for rice» and its principal producing Countries. 


Table A-2 tabulates the Countries and crops candidates for survey as 
seen from the UN/FAO outlook. 





APPENDIX B 


CURRENT CONVENTIONAL METHODS OF AGRICULTURAL CROP 
MEASUREMENT AND FORECASTING 

The measurement and forecasting of agricultural production are 
activities of major importance in virtually all countries of the world, 
as indicated in Figure B-1 [13] • 

The most sophisticated forecast system is the United States De- - 
partment of Agriculture (USDA's) Economic and Commodity Reporting 
Service (ERCS) Since much information is available on this system, 
it can conveniently be used to exemplify conventional procedures. 

The key functions of the ERCS data gathering. Information develop- 
ment and dissemination cycle are shown schematically in Figure B-2 [1]. 

Shown in Figure B-3 are SRS's three principal sources of data: 
a) voluntary reports from individual farmers (9.5 million question- 
naires mailed yearly, approximately 3 million responses); b) objective 
yield measurements by specialists; c) direct sampling of selected 
farms by specialists (enumerators) , through farmer interviews supported 
by map and aerial survey Information. 

A typical reporting form periodically filled out by individual 
farmers is shown in Figure B-4. Aerial photography at scale pf' order 
1:20,000 is used by USDA's enumerators to precisely mensurate a 
limited number (approximately 17,000, covering 0.6% of total farm 
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area) of agricultural producing units » The fanner's Information Is 
validated by planlmeterlng the farm's area, accurately Identified on 
the photograph. 

Yield prediction (bushels per acre or quintals per hectare) is 
performed by: a) the farmer, who assesses the "condition" of the 

crop (Its status with respect to what it was at the same time last 
year): and b) specialized personnel who perform "objective measure- 
ments" of various indicators of plant development on selected test 
plots; e.g., plant density, number of ears, number of splkelets, and 
so forth. 

These reports are Integrated by the central Crop Reporting Board, 
which issues monthly forecasts beginning 4 to 6 months prior to har- 
vest, depending upon the crop. A final estimate of production for 
the current year is given shortly after the harvest. 

USDA corrects their estimates over the following two years, based 
on additional information which becomes available on the amounts of grain 
products actually processed and sold by food Industries. As evidenced 
by Figure B-5, the corrections are quite small when compared with the 
variances of selected other countries. Figure B-5 Indicates that the 
fluctuations .In typical developing Nations reflect a considerably 
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higher degree of uncertainty than exists in developed Nations. 

USOA's method of forecasting production Is based upon the product 
of two factors; 1) the acreage planted; and 2) the estimated yield 
(bushels per acre or quintals per hectare). The yield estimate Is 
based upon a regression formula of the type: 
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FIGURE B-5 

VARIATION OF FINAL ESTIMATE 







B~6 


Y - aXj + bX2 + cXg + dX^ 

where: 

Xj^ is the "condition" of the' crop at time of estltiiace^ 
reported by the farmer or crop surveyor; 

X2 is the precipitation, in inches, which occurred during 
Che previous two months; 

is the precipitation, in inches, which is predicted will 
fall in the following two months; 

X^ is the historical growth of yield: for wheat, for example, 

approximately 3% per annum. 

a, b, c, d are coefficients whose value is derived from re- 
gressing historical records over approximately the pre- 
vious fifteen years. 

Implicit in this formula is the hypothesis that the methods of 
cultivation: a) fertilizer input; b) weeding effectiveness; c) den- 

sity of planting, and so forth, are essentially constant from year to 
year, and will be practiced by the individual farmer with the same dili- 
gence and with the best technology known at present. Also implicit is 
the assumption of a "normal" year as regards crop disease. The occur- 
rence of diseases such as wheat rust or corn blight, or unusual events 
such as floods, is factored in by the farmers and specialized survey 
personnel, who amend their estimate, upon detection, in the next monthly 
reporting period. 

Note that in the formulation above, the only corrector of the 
estimates is precipitation: rainfall is assumed to be the principal 
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driving function for yield* While this assumption has been found 
adequate for the U,S., it may not be optimal for Nations subject to 
different climatologies* France, for example, employs the additional 
variable of integrated temperature (degree-days) for their yield 
model. 

A measure of the accuracy of forecasts referenced to the "final" 
yearly estimate for the last several years is Indicated in Figure B-6. 

Figure B-7 depicts the average historical trend of forecast accuracies. 
Figure B-8 shows the average accuracy experienced over the last 15 
years for winter wheat, the major variety of U.S. wheat. 

The relative contribution to the forecast of the acreage and yield 
components is shown in Figure B-9 , which compares the fluctuations in 
U.S. wheat production with the variations in acreage and yield. Note 
that acreage fluctuations account for almost 50% of the total production 
variance , 

The USDA system is substantially followed, with varying degrees of 
sophistication, by the developed Nations. In developing Nations, simpler 
systems are employed. Table B-1, based upon FAO data, categorizes the 
crop reporting systems employed by the majority of the world's nations. 

Note that the complete system composed of farmer reports, plus objective 
yield measurement plus enumeration by specialists is used in general in 
those countries which possess the highest levels of living and literacy 
rate. As the GNP per capita diminishes, crop measurement systems tend 
to eliminate farmer reports and rely primarily upon sample measurements 
by specialists, The frequency and sophistication of the sampling procedures 
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COUNTRES INVOLVED IN CROP SURVEY ANALYSIS 
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also tend to decrease with decreasing GNP per capita. Table I >2 synthesis 
the world's current crop rt. siting systems Into three fundamental pategorles: 
Developing » Intermediate, Advanced. The agricultural data>-gatherlng network 
of the US Is shown in Figure B-10, that of m typical intermediate county 
in Flginre B-11. 

The basis for potential improvements in acreage measurement procedures 
lies in the fact that present methods, for budgetary reasons, only measure 
limited samples. Increased sampling, or, in the limit, wall-to-wall 
measurements, would improve accuracy; performance of the current sampling 
task without laborious field surveys would reduce costs. 

As regards yield estimates, the benefits of Increased sampling still 
apply but with significant differences. Expected yield is based to a 
significant degree upon the assessment of present crop condition. To the 
extent that this quantity can be measured more accurately by increasing 
the number of samples, yield forecasts will be improved. However, yield 
also depends upon other factors more difficult to estimate: principal 

among these is weather. Therefore, until the state of technology will al- 
low accurate prediction of weather, an irreducible uncertainty will remain 
in the forecast. Better sampling methods can, however, materially assist 
in reducing the errors of estimate. 






TABLE B-2 



THE THREE BASIC LEVELS OF CROP FORKASTING SYSTEMS 


Developing 

Inteisediate 

Advanced 


Typicals Ooainican Republic 

Typical: Italy 

Typical; United States of America 

Regional 
Adninistra- 
tive Struc- 
ture 

1599 sections grouped into 69 
onnunes 

7851 Covnmes grouped into 
91 Provinces 

1700 enwerated areas grouped into 44 States 

Crops Surveyed 

HaizOg riccg beanSg potatoesj 
onions g garllCg peanuts g co- 
conutSg oranges, bananas, cocoa, 
sugar cane, plantains, pineapple, 
coffee (in podl, avocado^ pears, 
cotton, tobacco 

Wheat, rye, barley, oats, 
maize, rice, sugar beets, 
potatoes, peas, beans, vine- 
yards, fruits, olives, lin- 
seed rapeseed, v^etables, 
tobacco, fibers 

Grains, fodder crops, tuber and root crops, 
sugar crops, pulses, oilseeds, hay and grass 
seeds, vegetable seeds, fruits, nuts, veget- 
ables, tobacco, fibers 

Methods of 
Data Gather- 
ing 

• Interview of producers by en- 
umerators 

0 Crop area from personal 
judgement, supplemented by 
cadastral survey 
• Cro:» yields from local in- 
qulry 

• Direct Inquiry to farmer respondents 

• Enumerators with aerial photos 

• Objective yield measurement by specialists 

Saapling 

Structure 

None 

[ Simple sampling procedure 

i 

Multi-frame stratified sampling procedures 

Organization 

Municipal Statistical Board in 
each commune 

a Data from commune col- 
lected by local correspon- 
dent assisted by ptcnrinciai 
agricultural inspector 
• Central Institute of St- 
atistics issues technical 
directives and publishes 
results 

1 • 5R5 HQ staff supported by 44 State offices 
' comprising 9 crop regions, 

• Refer to Figure 9 

Frequency 
of Crop 
porting 

Every 3 months 

Two crop reports per year 
First estimate at planting 
time; second estimate at 
harvest 

IhJltiple crop reports per year* Intentions 
to plant -yearly per crop. Acreage, crop con- 
condition, production forecast - monthly for 
3-6 months* Final production and yield-yearly. 












FIGURE B-11 

TYPICAL FOREIGN COUNTRY CROP INVENTORY & FORECAST INFORMATION NETWORK (ITALY) 


Acreage : 

Statistical Sampling 
(5,300 segments, 
stratified, 750 
enumerators) 



Provincial 


Inspectoi dtcs 


of Agriculture 


(92) 




Publlc-at-large 
National Farmers Association 
National Cattle Growers Co-op 
Chambers of Commerce 



Dealers and Traders 
AIMA (Seml-Gov*t, Purchasing 
Agency) 

ISVET (Private Forecasting 
Group) 

Private Traders 


Government Agencies 
Ministry Agriculture & 
Forestry 


L Yield: Subjective 

plus Farmer 
Merles (same per- 
sonnel) 


— Mall/Courler /Publications 
Conventional electrical 


Foreign Agencies 
USDA Attache' 
Other Attache's 


National 

Weather 

Service 


International Agencies 
FAO 
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DETERMINATION OF THE ACCURACIES OF AGRICULTURAL ACREAGE 
ESTIMATION ACHIEVED BY ORBITING DATA SYSTEMS 

(LANDSAT) 
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APPENDIX C 

DETERMINATION OF THE ACCURACIES OF AGRICULTURAL ACREAGE ESTIMATION 
ACHIEVED BY ORBITING DATA SYSTEMS 
(LANDSAT) 

X . 0 Accmracles Achieved by ERTS Investigators 

The ERTS crop classifications represent an Intensive series performed 
during the early phases of the LANDSAT program from its launch In July 1972 
until 1974. The spacecraft was at that time named ERTS, later changed to 
LANDSAT. 

The available principal Investigations reporting numerical results 
were used in the analysis [16] [17] [18] [19] [20] [21] [22] [23] [24]. 

Table C-l synopsizes 'the investigations analyzed. Nine ERTS investi- 
gators performed 224 distinct tests in ten locations, for seven major crop 
types. Thirty-two of these experiments Included sufficiently complete des- 
cription of the corresponding ground truth for complete comparison. Forty- 
-three provided aggregate results against various estimates of the ground 
truth e.g. SRS estimates. State Yearbook data etc. The balance Indicated 
only the percent of pixels correctly identified. 

The cumulative curve of the results, presented in Figure C-IA, graphs 
the accuracies achieved versus the percentage of the tests achieving them. 

It presents the results for both the inventory and the mapping, or land use. 


TAeU C-1 




SUMMARY OF ERTS-1 ACREAGE MENSURATION RESULTS 






AVERAGE 

METHOD OF 

RANGE 

INVESTIGATOR 

mi 

DATE 

MAJOR CROPS 

FIELD SIZE 

COMPUTATION 

OF ERROR 


Hill Co., 
Montana 
Holt Co. » 

8/7/7Z 

Wheat 

42,8 Ha. 


0.1-35S 

Ra Bryan Erb« 

Nebraska 

7/29/72 

Com 

366.3 Ha. 

Planimeter 

NASA (87) 

Imperial 
Co, ^ Ca 1 , 
Butte Co . 1 

11/6/72 

Asparagus 

45.9 Ha. 




Cal. 

9/19/72 

Rice 

255,2 Ka. 



0*S. S^monett 

Williams 


Spring Wheat 


Pixel Count 

0-2Sf 

Co., N.D. 

6/5/73 

Durum 

15. B Ha. 

(Proportion 



Oats 

Barley 


to Maximum) 






/^Vegetables 

No-64. 8 Ha. 



Cratg L. 

Hidalgo 

1/21/73 

J Citrus 

Pixel Count 

2-3SOL 

Ulegandi AR5 

Co., leK. 

5/27/73 

( Cotton 

So- 6.9 Ha. 

(31} 



1 Sorghum 

Cen-6.9 Ha. 





V^Rangel and 

Tot-12.2 Ha, 






rWheat 





Lubbock 

6/18/73 

\ Cotton 




ttarlon F. 

Co., Tex. 

7/24/73 

v-Sorghum 


Pixel Count 

1-35% 

Baumgardner, 

Greeley 


/Wheat 

NA 

(Proportion 

LARS (5) 

Co.f Kan. 

6/19/73 

V Pasture 


to HaxImuTi) 


Humboldt 


/Com 





Co., Iowa 

B/26/73 

V Soybeans 


Bias Corrected 


Oavtd Aa 

3 Counties 

8/9/72 

Com 

NA 

& Uncorrected 

11-65% 

Landgrebe, 

LARS ns) 

In 111. 


Soybeans 


ERTS results de- 
termined on 






basis of equal 







A non-equal 
weights. 



12 Counties 

8/21/73 

Corn 

NA 

Pixel Count 



In Ind. 


Soybeans 



Humberto C. 

Pampus of 

10/5/72 

Wheat 

NA 

Pixel Count 

1-55X 

Caruti . 

Argentina 


Barley 


Argentina (6) 



Rye 







Corn 




J.Pa Hahl- 

Ames, Iowa 

8/72 

Com 

NA 

Planimeter 

7-40% 

&tede, Iowa 

Doon, Iowa 


Alfalfa 



State (10) 



Soybeans 




Ralph Shay, 

Indiana 


Corn 

6-17 Ha. 

Pixel Count 

0.5-21% 

NASA (14) 

1144 nots 


Soybeans 


(Pro, to Max) 

Gene A. 

San Joa- 

7/2(772 

Range 

17 Ha. 

Pixel Count 

0-380% 

Thorley, 

quin Co., 


Vegetables 



UnIVa of 

Calif. 


Fruit 


(Proportion 


Calif. (ZZ) 





to Maximun) 



Numbers in parentheses Indicate the nunber of tests performed b^ each itrvesti gator. 



FIGURE C~1 

DISTRIBUTION OF CLASSIFICATION ACCURACIES FOR ERTS INVESTIGATIONS ON 

ACREAGE MENSURATION 


B. 





A. 


CLASSIFICATION ACCURACY OF ERTS INVESTIGATIONS 



NUMBER OF DATA 
POINTS: 224 


CLASSIFICATION ACCURACY OF ERTS INVESTIGATIONS WITH PROPORTION 

OF TARGET CROP >10Z 



NUMBER OF OA1A 
POINTS: 64 
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mode. Inventory mode means the determination of the proportion of a given 
crop within the area under examination« without however specifying the 
crop's geographic location. It is the mode universally used for estimating 
agricultural acreage. Mapping mode» as the word indicates, connotes the 
identification not only of the species sought, but also of their geometric 
boundaries and geographic locations. This mode is generally not employed 
in conventional agricultural acreage estimation, due its higher cost with 
respect to the inventory mode. 

Mote that the accuracies achieved in the mapping mode are lower than 
those competing to the inventory mode; this is quantitatively consistent 
with the fact that the mapping mode requires more Information than pro- 
portion estimation: conversely, for equal Informatlion content, its level 

of performance should be lower. 

The accuracy of multlspectral classification worsens when the propor- 
tion of target crops to the total scene is small ; a subset was extracted 
from the total statistical population, based on the criterion of 10% or 
larger proportion of target crops. This yielded 64 distinct data points: 
the test results are presented in Figure C-IB. Note the higher accuracies 
achieved. 

* 

A breakdown of the mode - whether inventory or land use - employed 
for each experiment by each ERTS investigator is given in Table C-2. 
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TABLE C-2 

BREAKDOWN OF THE NO. OF EXPERIMENT VARIANTS FOR EACH ERTS 
INVESTIGATOR BT MODE 

LANO USE 
PLUS 


INVESTIGATOR 

LMO USE 

INVENTORY 

INVENTORY 

TOTAL 

mmmt 

51 


6 

86 

m 

. 16 

0 

13 

2$ 

fiARirri 

1 

B 

0 

9 

UltGAMD 

4 

e 

0 

TO 

mm 

76 

0 

0 

76 

EARTHSAT 

0 

0 

1 

1 

KAHtST£0£ 

0 

0 

2 

2 

BAUMGARDNER 

1 

0 

0 

1 

THORLEV 

0 

0 

10 

TO 

TOTAL 

149 

43 

32 

224 


One-hundred slKty-nlne experiments reported the Individual field 
areas. The distribution of the field sizes shown in Figure C-2 indi- 
cates the predominance of fields below 50 hectares and above 10,000 
hectares. 

The sensitivity of the classification to area Is Indicated in 
Figure C-3 for each of the area classes given in Figure C-2. The general 
trend is that the mean error within each class is Inversely related' to 
the area. Note however that there is a wide dispersion of errors for 
each area class: this is indicative of the wide dispersion in the frac- 

tion of the total scene comprised by the crops investigated. Figure C-4 
summarizes the results for the 169 investigations. 
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FIGURE C-3 

SENSITIVITY OP EFTS CLASSIFICATION TO AREA - INVENTORY MODE 

169 DATA POINTS 


0-50 ho. 



'ia-20 *20-50 ' >50* 


<1 1-5 5^0 



100 - 250 ho. 




500- 1,000 ha 




% ERROR 


5^000- 10,000 ha. 


c I 



T 


1 



1-5 ’5- >0 *10-20 '20-50 >50 


% ERROR 


1-5 5-tO ' tO-20 20-50 ^50 




The effect of small crop acreages was eliminated by reducing the 
set of Investigations analyzed from 169 to 64 as Indicate! Ln Table C-3. 
Figures C-5 and C-6 present the distributions of Investigations as a 
function of percent error for the 64 investigations with acreage frac- 
tion of crop greater than lOJt!. 

TAOIC c-3 

AGRICULTURAL CROP CLASSIFICATION FROM CRTS INVESTIGATIONS REDUCED TO 
PROVIDE ACREAGE FRACTIW > TO! 

NUHOER OF INVESTIGATIONS SELECTED FOR THIS SAMPLE: 32 WITH 64 CROP DATA 

POINTS 

NUMOER OF INVESTIGATORS PERFORMING THE 32 INVESTIGATIONS ANALYZED: S 

RANGE OF AREAS CLASSIFIED: 400 TO 100.000 HA: MEDIAN 1700 

GEOGRAPHIC LOCATIONS; 9 STATES: (IL, IN. HT. lA. CA. NE, GA. TX. NO); 

APPROXIMATELY 20 COUNTIES 

PRINCIPAL CROPS SELECTED: WHEAT. SHALL CRAINS, CORN, SOYBEANS. BARLEY, 

PASTURE , ALFALFA 

AVERAGE NUMBER OF DISTINCT CROPS CLASSIFIED PER INVESTIGATION: 3 TO 4 

NUMBER OF DISTINCT CROPS SELECTED FOR COMPARISON: AVERAGE 2 PER INVES- 
TIGATION. CRITERION; ONLY CROPS WHOSE ACREAGE RATIO TO TOTAL CROP 

ACREAGE EXCEEDED 10! 

THESE YIELDED 64 DATA POINTS ON PRINCIPAL CROPS. 

The classification results are presented In Table C-4 for both 
the 224 and 64 investigation data sets. It shows that the one sigma 
accuracy of classification In the ^nventoAt/ mode, for the total popula- 
tion of 224 data points is 74%: this Increases to 86% for the subset 

of 64 data points featuring 10% or greater proportions. 





ALL CROPS >10:S OF SCENE 



FIGURE C-5 
LAND USE MODE 
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FIGURE C-6 
INVENTORY MODE 
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TAeiE C>4 

ACCURACY OF CLASSIFICATION OF CROPS ACHIEVEO 
BY CRTS INVESTl^TORS 


UNIVERSE OF 224 CLASSIFICATION TESTS ON 7 CROPS, 10 GEOGRAPHIC LOCATIONS, 
ONE SIGMA CONFIDENCE 

PROPORTION 

ESTIMATION WPPING 

ACCURACY X 74 G3 

ERROR X 26 37 

64 aASSIFlCATION TESTS EXTRACTED FROM THE ABOVE UNIVERSE, WHEREIN THE 
CROP-TO-TOTAL-AREA PROPORTION HAS GREATER THAN lOX. 

ACCURACY I 86 76 

ERROR t 14 24 


2,0 Accuracies Achieved In Acreage Estimation after 1974 

ERTS investigations ceased after 197A. To assess whether the early 
performance hna improved in tbe approximately four years- elapsed since, a 
number of subsequent investigations, performed by reputable organizations, 
were analyzed. Table C~5 is a summary of the key Landsat investigations 
dealing with crop classifications reported since 1974. 

Figure C-7 presents the results of the investigations in the Inventory 
and land use mode. The one sigma accuracy was 86.3% for the inventory mode 
and 73.4% for the land use mode respectively. The overall results are also 
shown in Table C-6, 
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TABLE C-5 



. 


suwwir or kev iahosat classification investigations since igta 



0£StRIPT10N 

LOCATION 

NUWEIt OF 
CLASSES 

AVERAGE 

INVENTORV 

ERROR 

AVERAGE 
LAND USE 
ERROR 

DATE 

PROCEDURE 

Are« cstittatfon of 
crops by dfglU] 
analysis of 
Landsat Data. LARS, 
Purdue University 

Kansas 

4.613.000 Ha. 

2 

(wheat, other) 

1.51 

compared with 
SRS estimate 

— 

1977 

Computer classifi- 
cation of wheat. 

Klgh-altitude v$. 
Landsat Inagery for 
digital crop Iden- 
tification (LACIE) 

Kern County, 
California 
11,839 acres 

10 

(cotton, grapes, 
melons, toma- 
toes, sugar- 
beets. wheat, 
fallow, oranges, 
natural vegeta- 
tion. other) 

24X 

22 % 

1977 

Computer crop 
classification 

Pattern Recognition 
of soils a crops 
from space » USDA 

Imperial 
Valley. 
California 
(S3 fields) 

5 

(alfalfa, bar- 
ley, sugarbeets, 
bare soil, salt 
flats) 


2B% 

197S 

Computer crop 
classification 

land use inventory 
for Alto Lazio, 
ERIH 

Alto Lazio, 
Italy 

(53.264 ha.) 

'S 

(Piscolo, Semina- 
t'fvl, Foreste, 
Orchards. Urban) 

12,3% 


1977 

Manual crop 
class from 
Landsat Images 

Landsat spectral 
signatures: 
studies with soil 
associations and 
vegetation. South 
Dakota Agricultural 
experiment station 

Brookings Co. 
South Dakota 
12.950 ha. 

4 

(com, small 
grain, grass, 
water) 

6.U 

6.3t 

1978 

Computer clas- 
sification of 
landsat Imagery 

Cost, accuracy, 
and consistency 
•comparisons of 
land use maps 
made from high- 
altitude air- 
craft photography 
and ERT5 imagery. 

Central 

Atlantic 

6 

(urban, water, 
agricultural . 
forest, non- 
forested wetland 

12.U 

30.5t 

1975 

Manual classifi- 
cation from 
Landsat imagery 

Landsat urban 
area delineations. 
Incralab Project 
75-3 

Prince 

Georges Co.^ 
MD 

64,100 ha. 

2 

(residential, 

other) 

27 % 

27t 

1977 

Computer classi- 
fication 

Global agricultural 
productivity evalua- 
tion from Landsat 
data, J. Schubert, 
Gregory Geoscience 
Ltd. 

4 Canadian 
test sites 
and acres 
analyzed 
1) Swiftcur- 
rent 

(1644 acres) 
(7204 acres) 

2 

(wheat, other) 

0 % 

0.6S 


1976 

July 

1973 
July 

1975 
July 

1976 

July 

1974 
July 
1976 

Aug. 

1976 

Manual classifi- 
cation from en- 
hanced Landsat 
imagery 


(6960 acres) 


3.4X 

- 



2) Raymond 
(2175 acres) 

2 

(wheat, other) 

3.8% 

- 



(4611 acres) 


-4.7X 

- 



3) Helford 
(4206 acres) 


-I3.lt 

/ 

- 




-•a 

1 





TABLE C-5 (cont*<0 


SUHHARY OF KEV LAKOSAT CLASSlflCATlON INVESTIGATIONS SINCE 1974 


DESCRIPTION 

LOCATJOPi 

NUWER OF 
CLASSES 

AVERAGE 

INVENTORY 

ERROR 

AVERAGE 
LAND USE 
ERROR 

DATE 

PROCEDURE 


4) Stoney 
Mountain 
(3439 
acres) 

2 

(wheat, other) 

7*2X 

— 

Aug. 

1974 


Crop Identification 
tethnoiogy assumed 
for remote sensing 
(CITARSh B.d. 

Oavis, A.H. Eeiueson, 
NASA, JSC, Houston, 
TeKas 

Texas 

3 

(com, soy- 
beans, other) 

13. 2X 


1975 

Computer classic 
ficatlon 

Land use classifi- 
cation for 6 
Rocky Mountain 
states - using 
Landsat multispec- 
tral - 
poral data 

Fox Creek, 
Colorado 
(1268 
acres) 

n 

(Rabbitbrush, 
grass/rabbtt- 
brush, dense 
shrub, mahogany, 
meadow, wet pas- 
ture, doug fir, 
Pinyon-Juniper, 
Ponderosa Pine, 
Aspen, cotton- 
wood 


25. 6S 

1977 

Computer classi 
fication 


Fox Creek, 
Colorado 
("■'68 
acres) 

5 

(shrub land, 
non- irrigated 
grassland, ir- 
rigated grass- 
land, coni- 
ferous forest, 
deciduous for- 
est) 


17. 6X 

• 



2 

(Rangeland, 


}OnS% 



FIGURE C-7 






■^^LANDSAT POST 74 CROP CLASSIFICATION RESULTS 


TOTAL NO. IN SAMPLE - 18 



1 SIGMA INVENTORY ERROR :I3.7% 

1 SIGMA LAND USE ERROR Z26.6% 

MEDIAN ABSOLUTE INVENTORY ERROR Z 8 % 
MEDIAN ABSOLUTE LAND USE ERROR Z23% 
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. T*BI,£ C»6 

ACCtmCY OF CLASSIFlCATtOH OF CROPS ACHIEVED FRM 197< TO 1978 

Universe of 16 classification cnUeavoirs on ID cropSt 13 geographic 
locations, one slgna confidence. 


PROPORTION 

ESTIMATION HAPPING 


Accuracy i 

67% 

73% 

Error X 

m 

27% 


<JKlGli\AL PAGE LJ 
OF POOR QUALITYi 


3.0 Description of the Ex p eriments In ERTS Crop Classification 

Tables C-7 through C-14 present the highlights of each ERTS experiment: 
the crops analyzed, their locations, the dates of the ERTS observations, 
the method of analysis employed. 

Figures C-8 through C-16 depict the scattergrams of the differences 
between measured values and ground truth. The scattergrams present for 
each group of investigations, Che error experienced in measuring each crop" 
area versus the actual areas. .The solid downward sloping straight lines 
represent the error of mensuration e^% expected from the ERTS pixel size, 
on Che assumption of perfect discrimination. They obey-' 'the relationship; 


/~A~ 

* 

r Is the average dimension . of the ERTS pixel Z70 meters 
A Is the area mensurated In hectares 




Where: 
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TABLE C-7 


HIGHLIGHTS OF ERTS EXPERIMENT 




^ ■ 

;:i j 


Analysis, R. Bryan Erb, 

NASA, 

August 1974. 



Scene Content: 





Hill County 

Summer Fallow & Sod 

55% 

Imperial County 

Bare Soil 

40% 


Spring Wheat 

24% 


Asparagus 

34% 


Winter Wheat 

19% 


Alfalfa 

15% 


Barley 

2% 


Carrots 

5% 





Melons 

4% 





Light Veg. 

2% 

Holt County 

Field Com 

67% 

Butte County 

Rice 

52% 


Grass & Pasture 

16% 


Fallow 

41% 


Sunflowers 

4% 


Com 

6% 


Popcorn & Field Com 4% 


Alfalfa 

1% 


Popcorn 

4% 


Sugar beets. 5% 


Field Com & Grain 






Sorghum 

3% 





Disked Fallow 

2% 




Time of Measurement: 





Hill County 8/7/72 


Imperial County 

11/6/72 



Spring Wheat - Ripening 
Winter Wheat - Harvesting 
Barley - Harveatipg 

Holt County 7/29/72 

Field Com - Tasseling 
Grass & Pasture (Native 
Hay) “ Harvesting 
Sunflowers - Heading 
Grain Sorghum - Heading 


Butte County 


Alfalfa - Growing 


9/19/72 

Rice - Maturing & 
Harvesting 
Com - 100% Ground 
Cover 

Sugar beets - Harvesting 


Imperiar County - 45 « 9 Ha. 
Butte County - 225.2 Ha. 


Average Field Size ; 

Hill County - 42.8 Ha. 

Holt County - 366.3 Ha. 

Method of Computation ; 

Mechanical/electronic planlmeter to manually trace field Jiorders. 

Additive color-enhanced ERTS-1 Imagery of about 1:130 >000 scale used as base 
photography for measurements. 

Universal Transverse Mercator coordinate system with grid prepared for en- 
hanced ERTS-1 imagery of about 1:130,000 scale to locate fields. 












TABLE C>8 




HIGHLIGHTS OF ERTS EXPERIMENT 


Source : Reflectance of Vegetation, Soil, and Water , Craig L, Wlegand et al*. 

Agricultural Research Service t USDA« November 1974. 

Scene Content : 

Urban 

Agricultural 

Vegetables 

Citrus 

Cotton & Sorghum 
Idle Cropland 
Dry Debris 
Rangeland 
Grass 
Shrub 

Non-agricultural 

Water 

Time and Place of Measurement ; 

January 21 and May 27 f 1973, Hidalgo County, Texas. 

Average Field Size : 

Total County - 12.2 Ha. 

Northern Region Only (Range and Pasture) •64.8 Ha. 

Central Region Only (General Faming & Citrus Growing) - 6.9 Ha. 

Southern Region Only (Winter Vegetables) --8,9 Ha. 

Method of Computation ; 

Ground truth information was obtained on several Items pertaining to crops, 
crop histories, and faming practices by interviewing approximately 400 
famers. Typical ground truth included crop species, stage of maturity, 
row spacing, plant height, percent ground cover, row direction, Munsell 
color of soil, recent cultural practices, and other qualitative information. 

A classification was determined for every pixel in Hidalgo County (849,000 
pixels in January and 948,000 pixels In May), and a comparison was made, 
using Student's t-test, between actual and computer acreage estimates of 
Hidalgo County for the categories of vegetable, citrus, cotton and sorghum, 
idle cropland, grass, and mixed shrub. 
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TABLE C-9 


HIGHLIGHTS OF ERTS EXPERIMENT 


Source ; Evaluation and Comparison of ERTS Measurementa of Major Crops 
and Soli Aasoclntlons for Selected Test Sites In the Central 
United States . Marlon F. Baumgardner and Staff, LARS, Purdue 
University, March 1974. 

Scene Content ; 

Lubbock Co., Texas - Wheat " harvesting 

Cotton 

Grain Sorghum 
Hater 

Greeley Co., Kan. Ha..'d Red Hlnter Wheat - heading 

Fallow (Bare Soil) 

Permanent Pasture 

Humboldt Co., Iowa - Com 

Soybeans 
Other Row Crops 

Time of Measurement ; 

Lubbock Co., Texas - June 18 and July 24, 1973 

Greeley Co., Kan. - June 19, 1973 

Humboldt Co., Iowa - August 26, 1973 

Method of Computation ; 

** 

Lubbock Co., Texas -* Acreage estimates were made from percentage of points 

in each class. The total area In the county was mul- 
tiplied by these percentages to obtain an estimate of 
area In each cover type. 

Greeley Co., Kan. - Area of wheat was calculated from classification results. 

Humboldt Co., Iowa..- Area of crops was calculated from classification results. 





FIGURE C-11 


EXPERIMENTAL ERRORS IN GEOMETRIC MEMSURATION FROM U^SAT IMAGERY ■ 
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TABLE C-10 


HIGHLIGHTS OF ERTS EXPERIMENT 


Source ; A Study of the Utilization of ERTS-1 Pata from the Wabash River 
Basin , David A. Landgrebe and Staff, LARS, Purdue University, 
August 1974. 


Scene Content: 


Illinois - 


Com 40% 
Soybeans 18% 

Grain Sorghum 

Alfalfa 

Hay 

Pasture 

Small Grain Stubble 

Woods 

Urban 


42% 


Indiana - Com 32% 

Soybeans 25% 

Hay 

Pasture 

Small Grain Stubble 

Woods 

Urban 


43% 


Time and Place of Measurement ; - 

Illinois - 8/9/72 
DeKalb County 
Ogle County 
Lee County 

(Analysis performed on total Information 
from these three counties.) 


Indiana - 8/21/73 
Lake County Newton Cbunty 
Porter CoOnty White County 
LaPorte County Benton County 
Starke County Warren County 
Pulaski County Fountain County 
Casper County Tippecanoe County 
(Analysis performed on total in- 
formation from these twelve 
counties.) SEE ATTACHE]^ MAP. 


Method of Computation : 


Illinois - Bias corrected and uncorrected ERTS- results were determined, each 
on the basis of equal weights and non-equal weights, resulting 
in four acreage estimates made from ERTS data. The method used 
for unblasing classification results Involves multiplying the 
county classification results by the inverse^ of the test field 
classification performance matrix as' follows^ 


,-l 


Where: 


CP 


A B 1 X n vector of crop acreages 

^ classification vector with n crops or classes 
- inverse of n x n classification performance 
matrix 


Indiana - Acreage measurements based on pixel count classification. 
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TABLE C-11 


HIGHLIGHTS OF ERTS EXPERIMENT 


Source: Identification and Quantification of Agriculture Crops, Deter~ 

mlnatlon of Dynamics of Land Use ond Agricultural Census , Eng. 
Humberto C. Garutl, Secretarla De Agriculture Y Ganaderla, 
Argentina, February 1974. 

Scene Content ; 

Wheat - eared (flowering - grained) 

Barley 

Rye 

Com - geminated 

Grain Sorghum " ready to plant 

Forage Sorghum - ready to plant 

Soya •> planted 

Mljo - planted 

Sunflower •> planted 

Alfalfa - shoo ted 

Mixed Pasture ^ shooted 

Natural Pasture - shooted 

Winter Cereal 

Time and Place of Measurement ; 

Pampus Region of Argentina - 10/5/72 
Method of Computation ; 

1. Manual Interpretation: Acreage estimates were obtained from 1:1,000,000 

amplified transparencies of false color images after classification. 

2. Pixel Count: Ground truth available for only one scene; since the 3 

scenes used were successive frames from one ERTS pass, ground truth used 
as training data from first scene, then processed results from first 
scene used as initial training data for next scene, etc. (10% overlap 
between scenes. 




- 





Conclusions : 


The best method for Investigating large cultivated areas Is automatic 
discrimination* which assures speed and consideration o£ large areas. 

A 96% safety factor was obtained in the general preliminary classifica- 
tion. 


TABLE C-12 


HIGHLIGHTS OF ERTS EXPERIMENT 

" II. ■ ^ m y, 1,^. , 


Source ; Remote Sensing In Iowa Agriculture; Identification and Clas- 
sification of Iowa*e Crops, Soils, and Forestry Resources 
Using ERTS-1 and Complimentary Underflight Imagery , J.P, Mahl- 
stede et al.» Iowa State University, May 1974. 

Scene Content ; 

Com 

Alfalfa 

Soybeans 

Oats 

Pasture 

Water 

Forest 

Time and Place of Measurement ; 


Ames, Iowa *- August 1972 
Doon, Iowa - August 1972 


Method of Computation ; 
Optical Processing 
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TABLE C-13 

HIGHLIGHTS OF ERTS EXPERIMENT 

Source ; EOS Thematic Happer Technical Review . Ralph Shay et al 
April 30 - May 2, 1975. 

Scene Content ; 

Com 

Soybeans 

Time and Place of Measurement ; 

Midveat (Indiana and Illinois 
Average Field Size ; 

6-17 Ha. 

Method of Computation ; 


Pixel Count 
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TABLg C-14 


HIGHLIGHTS OF ERTS EXPERIMENT 




Source : Regional Anrlcultural Surveys Using ERTS-1 Data» Gene , 

et al*» Center for Remote Sensing Research, University 
nla. UN 640. 

Scene Content: 


Range 

Pasture 

Plowed Land 

Grapes 

Alfalfa 

Barley 

Sorghum 

Tomato 


Sugar Beets 

Walnut 

Fruit Trees 

Burned Areas 

Beans 

Harvested 

Cucumber 


~Tlme and Place of Measurement; 


San Joaquin County, California July 26, 1972 


Average Field Size ; 

Approximately 17 Ha. 

Method of Computation : 

Manual Interpretation 
Automatic Interpretation 


i 


L, Thorley 
of Califor- 


r r ” ~y;;; 

FIGURE C~16 

EXPERIMENTAL ERRORS IN GEOMETRIC MENSURATION FROM LANDSAT IMAGERY - THORLEY 
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APPENDIX D 

ANALYSIS OF AVAILABLE CROP SPECTRA 

It has been known for a long time that laboratory spectra of single 
leaves exhibit a considerable degree of qualitative similarity. Personal 
communication by Professor R.N. Colwell of the University of California, 
Berkely, In 1975* who obtained numerous laboratory spectrograms* indicated 
that essentially "all leaves look alike." An example of this qualitative 
likeness Is depicted In Figure D-1 which shows the qualitative similari- 
ties In the reflectances of single leafs in the spectral region from 
0.5 to 2.5 vm of various plants: Bean, Avocado, Sorghum, Pigweed, Com, 

Soybean, Orange, Peach, Cotton. 

1 , 

Leaf spectra however are not necessarily good indicators ' of the spectral 
reflectances of entire plants: especially since the available leaf spectra 

are taken "In vitro" at essentially normal incidences. One could speculate 
that specific patterns of leaves peculiar to each plant may emphasize pre- 
ferred directions of reflection, or possibly even enhance one spectral region 
over others. 

Analyses of spectra contained in the literature, see Table D-1, did 
qualitatively confirm the hypothesis of "similarity" of spectral shapes: 
however, the difficulties encountered In resolving discrepancies of cali- 
bration, stage of plant development. Illumination and viewing angles, 
prohibited the formulation of general relationships valid between diverse 

plants, even between plants belonging to the same species but at differing 

4 '- ■ 

stages of development, or subjected to different stresses. 
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TMU M 

nmcipUL »Ecnwi mta mscs owiiwie iwjg litctatwe 


KSIOUTION 

SPONSOR 

R6ENCV 

MEASUREWKT 

DATE 

LOCATION 

KTF 

MSR/NTF 

W73 - 5/74 

KTF - Shy St. Loult, NS 

IMS 

MASR 

7/68 

LARS • Vest Lifaytttt. IN 

mh 

U.S. AIR FORCE 

1972 

NTIS > Sprlngftjild, VA 

ERDL 

tt.S. ARNI 

1961 - 1962 

NTtS - Sprfngfteld, VA 

WSMO 

USM/MtS 

]0/66 • 9/67 

USDA ARS • Htilico, TX 

XRINOV 

USSR 

1936 « 1938 

NRC - OtttMt. Cituda 

IITRI* 

NASA/USAF 

ViHous 

HTIS • Springfield, VA 


*SwMry of M» from other toums. 


The resolution of this question had to await the advent of flne-pralned, 
calibrated spectra. 

Figures D-2 and D~3 present spectral radiance values produced by the 

Goddard Institute for Space Stqdles from a fl^ld program conducted In 

_ * 

Imperial Valley » California* Figure D-2 shows spectra of the crops al- 
falfa» melons and sorghtim at a 40 to 50% leaf cover. Figure D-3 shows 
the spectra of the crops wheat, alfalfa, sugar beets, sorghum, asparagus, 
each at an approximate leaf cover of 100%. These plots show that diverse 
crops at similar stages of development, l.e. at similar amounts of leaf 
cover, exhibit qualitatively similar spectral signatures. 




400 500 600 700 800 900 1000 

WAVELENGTH" (nancjmeters) 


400 500 600 700 800 900 1000 

WAVELENGTH (nanometers ) 


400 500 600 700 800 900 100 

WAVELENGTH (nanometers) 


MELONS 


SORGHUM 


ALFALFA 


SOURCE: Atlas of Selected Crop Spectra IiBperlal Valley, California 

S.G. Ungar* NASA Institute for Space Studies* 1977 



PIGURk D-'S 

SPECTRA OP VARIOUS CROPS AT 100% ICflP COVER IM IMPERIAI, VALUBY, CALIFORNIA 


mmm 


400 500 600 700 800 900 1000 

WAVEIEMGTH (nanometers) 
WHEAT 


400 500 600 700 800 900 1000 

WAVELENGTH (nanometers) 

SUGAR BEETS 


400 500 600 700 800 900 lOOC 

WAVEIENGTH (nanometers) 

SORGHUM 



400 500 600 700 800 900 1000 

WAVELENGTH (nanometers) 

ASPARAGUS 



500 600 700 800 900 1000 

WAVELENGTH (nanometers) 

AIFALFA 


SOURCE i Atlas of Selected Crop Spectra Imperial Valley, California 
S.G. Ungar, NASA Institute for Space Studies, 1977 
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A more quantitative determination of the spectral similarity between 
crops at like levels of leaf cover 1s obtained by examining the ratio 
between various spectral bands. For example » Table D-2 presents results 
of a correlation analysis between the spectral bands 0.7 to 0.8 ym and 0.8 
to 0.9 ym for wheat and alfalfa crops with greater than 602 leaf cover. 




TABLE 0*2 

1 CORJIELATIW AHM.YS1S BEHIEEN BANOS 0*7«0.6 im AND O.B^O.9 urn 1 


Number of Spectra Used; 149 

RESULTS: 

COEFFtClENT 
OF CORREUTION 

Sn. ERROR 
SLOPE OF SLOPE 

VHEAT 

0,93 

0.58 0.05 

ALFALFA 

0.97 

O.SB 0.02 

soil 

0.98 

1.14 0.04 

Source: 

CCOsyiloRS analysis of a selected subset of spectra 
provided by Dr. S. Unger of the Goddard iRstltute for 
Space Science. 


Each crop exhibits a "slope’* (ratio of amplitude of band 0.8~0.9 to band 

, t 

0.7*-0.8 ym) of 0.58 and a high degree of correlation between the two spec- 
tral bands. The closeness of the slope values indicates that these crops 
generate similarly shaped spectra: whereas the high correlations Indicate 

that the computed slope, hence the shape of the spectra Is relatively 
stable for each crop at leaf cover areas greater than 60%. A number of 
such correlations were performed by ECOsystems with results similar to 
those obtained by other workers. 

Table D-3 lists relative reflectances of irrigated wheat in the two 
spectral bands "'0.74-0. 76 ym and 0.78-0.80 ym and the ratio between the two 




T ^ ' 
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for the phenologlc stages from early jointing through dead ripe for Bushlandi 
Texas* [ 3 These two bands show significant amplitude changes throughout 

the growth period. They exhibit 1 sigma variations of approximately 20% 
in the band 0.74-0.76 pm and 26% in the baud 0.78-0.80 pm. The ratio be- 
tween the two bands remains relatively constant* however* through this 
growth period, exhibiting a 1 sigma variation of approximately 8%. From 
these results it can be concluded that while the wheat field exhibits sig- 
nificant amplitude changes in mean reflectance throughout the later stages 
of the crop's growth, the shape of the spectra remains relatively constant. 




TABLE 0>3 




RATIO OF RELATtVC REFLECTANCE 0,785 um/0.74S i« FOR IRRIGATED 


UKEAT FiaO NO. 9Q10UI. 6USHLAN0. 

Jl 

* 


REaECTANCE 

REaECTANCE 

RATIO 

STAGE 

DATE 

3/16 

0.74>0.76ua 

0.78*0. 80pm 

jii_ 

.338 

.387 

1.145 

Eorly Oolntlns 

,331 

.382 

1.154 



.331 

.370 

1.118 



.315 

,337 

1.070 


4/4 

.376 

.520 

1.383 

Late Jointing 

.256 

.299 

1.16B 



.375 

.531 

1.416 

• 


.393 

.546 

1.389 


4/21 

.361 

.494 

1.368 

Booting 

.394 

.513 

1.302 



.387 

.498 

1.267 



.366 

.457 

1.249 


5/19 

.304 

.382 

1.257 

Early Ripa 

.298 

.373 

1.252 



.316 

,379 

1.199 



.302 

.349 

1.156 


5/27 

.248 

.304 

1.226 

Late Ripe 

.343 

.364 

1.061 



.229 

.245 

1.070 



,227 

.287 

1.264 


6/4 

.231 

.269 

1.165 


.222 

:MV 

1.158 



.227 

MSO 

* 


.241 

.277 

1.149 


HEAK 

.309 

.378 

1.215 


0 

.060 

,097 

.103 


c,v. 

.195 

.256 

.084 




This relationship appears to hold for other crops as well. Figures D-4 
end D-5 present spectral plots of various crops throughout their growth for. 
Hand County* South Dakota [ ] and Bushland* Texas. The similar ity of 

shape* although with significant variations in amplitude from stage to stage 
are evident In these plots. 

Figure D-6 plots single-date "snapshots" of spectra of various ground 
covers taken at the Bushland, Texas site. Note the emergence of irrigated 
wheat from the soil cover spectrum on the earlier date: its further growth* 

and the emergence of dryland wheat* on the later date. On the later date, 
whereas the spectral shapes of the emerging crops are qualitatively similar* 
their average magnltudi^s begin to differ significantly. 

A striking example of the differences of amplitude displayed by crops 
at certain times during their growth cycle is shown in Figure D-7. Note 
the significant emergence of the crop rye over the others: its spectral 

shape* however* is qualitatively not very different from the spectral 

shapes of other crops. It should be noted that the April 4, 1974 curve for 

rye appeared upon investigation to be probably tainted by an imperfection 

in the experimental procedure possibly an instrumental gain change of as 

much as 1.5:1. Even correcting for this possible mishap does however not 
tnodlfy the general relationships. 

Figure D-8 plots the temporal progression of the mean reflectance of 
spring wheat* pasture* com* oats* alfalfa and barley in the 

spectral band 0.76-0.9 pm for Hand County* South Dakota. It shows that 

temporal variations in mean reflectance differ from crop to crop. 


BIOIRCCTIONnt RSFLCCrnNCC FACTOR I'/.) 



BIOIRCCTtONRL REFUCTPNCC FACTOR (^1 



REFLECTRNCE OF SPRING WHEAT 
AT DIFFERENT MfiTURITY STAGES 


LOCnnON* HAND COUNTY. SOUTH OfVCOTR 
SENSOR' FSS DATE' 1976 



URVELENGTH HliaOICTERS) 



FIGURE D-5 













FIGURE D-7 
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It can be argued at this point that the findings thus far exemplified 
are peculiar to the limited number of regions within which they have been 
gathered: and that the hypothesis of "similarity" of spectral shapes Is 

thereby not proven, but rather only Inferrable as a working hypothesis. 

There Is of course no positive answer to this argument, whose thorough 
experimental verification would require an Impracticable number of field 
data gathering programs. Heuristic arguments which point to the fact that 
the hypothesis of "similarity" Is valid worldwide are: 

1. Personal communication by Dr. N.G. Khar In, Director, Remote Sensing 
Institute, Academy of Sciences of the Uzbekhlstan Republic, USSR. 
Or. Kharln confirmed -In 1976 that from his findings In the USSR, 
spectra of trees and crops appear "very similar." His tabulations 
in [25] further confirm this fact. 

2. Spectra of rice gathered In hothouse environment by the French- 
Itallan cooperative program "Agreste" show not only that various 
species of rice present "similar" reflectance spectra, but that 
these are essentially "similar" to the available spectra of other 
U.S. crops [26]. 

The mean reflectances of crops display fluctuations. Figure D-9 
shows the fluctuations of mean reflectances measured among five spring 
wheat fields on July 28, 1976 In the Williams County, North Dakota LACIE 
Supersite. It also plots the results of 5 spectral scans of spring 
wheat field 334 located within this site; It shows that fluctuations oc- 
cur not only among fields but within single fields. The Figures show that 
such fluctuations occur In all regions of the spectrum from 0.4 - 2.4 ym. 


BIDIRECTIONAL REELECTANCE FACTOR (F.l 


FIGURE D-9 


VflRinBILITY IN REFLE27=.;JCE 
WITHIN fl SPRING WHEfiT FIELD 


VRRIRBILITY IN R5FLFCTRNCE 
RHONG SPRING WH'IRT FIELDS 


LOCATION* HILL inns COUNTY, NORTH DAKOTA 
SENSOR* FSS DATE* JULY 28. 1976 

FIELD* 334 


LOCATION* MILUAMS COUNTY. NORTH DAKOTA 
SENSOR* FSS DATE* JULY 28. 1976 


SCAN ONE 

sew T«0 

sew TWEE 

SOT FOUR 

SOT FIVE 
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Figure D-IO plots coefficients of variation as a function of spectral j 

wavelength for wheat in Bushland, Texas. It» and several other similar | 

.4 

tests, indicates that an average spectral variability of order 10% Is ’ 

1 

typical for reflectance fluctuations within a single field. Tables D-4 

t. 

through D-6 list coefficients of variation for various multl~field crop ^ 

classes for six dates In Hand County. South Dakota. They Indicate that 1 

I 

multi-field coefficients of variation are generally greater than 10%. 

The maximum multi-field reflectance fluctuation Is exhibited by oats with 
a coefficient of variation of approximately 40% on June 16. 1977 In the 
spectral band from 2.08 to 2.35 pm. Typical coefficients of variation 
for multi-field crop classes appear to fall In the range from 10% to 30%. 






TABLE D-4 


SPECTRAL 

BAND 

0.450 - 0 
0.520 - 0 
0.630 - 0 
0.760 - 0 
1.550 - 1 
2.080 - 2 




COEFFICIENTS OF VARIATION OF REFLECTANCE VERSUS BAND AND DATE: HAND COUNTY. S.D. 


SPRING WHEAT CV% WINTER WHEAT CV% 



4/21 

5/10 

6/01 

6/16 

7/07 

7/27 

4/21 

5/10 

6/01 

6/16 

7/07 

7/27 

.520 

13.5 

7.1 

13.6 

15.8 

15.0 

15.0 



20.7 

33.9 

19.3 

22.0 

.600 

15.5 

6.6 

9.4 

14.0 

14.2 

16.9 



14.6 

31.1 

17.2 

18.4 

.680 

13.9 

6.4 

9.9 

14.7 

16.2 

18.0 



14.4 

34.0 

15.6 

16.2 

.900 

13.6 

6.7 

8.5 

13.6 

12.0 

16.6 



13.5 

29.0 

7.6 

14.2 

.750 

6.1 

3.5 

9,2 

16.0 

10.7 

7.1 



19.7 

39.7 

15.0 

12.6 

.350 

7.3 

4.8 

17.5 

25,8 

19.4 

8.5 



34.1 

48.4 

20.5 

21.4 


AVERAGES OF 5 FIELDS, 74 TO 148 SPECTRA 
TAKEN AT EACH DATA 


AVERAGES OF 4 FIELDS, 50 TO 109 SPECTRA 
TAKEN AT EACH DATE 
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TABLE D-5 


SPECTRAL 

BAND 

0.450 - 0 
0.520 - 0 
0.630 - 0 
0.760 - 0 
1.550 - 1 
2.080 - 2 


COEFFICIENTS OF VARIATION OF REFLECTANCE VERS 


CORN CV% 



4/21 

5/10 

6/01 

6/16 

7/07. 

7/27 

.520 

15.5 

20.1 

8.6 

9.7 

5.2 

13.9 

.600 

16.6 

18.7 

8.1 

10.4 

5.0 

8.7 

.680 

15.1 

18.0 

8.4 

10.6 

6-0 

9.5 

.900 

12.9 

23.5 

8.4 

7.0 

1-8 

1.8 

.750 

6.8 

14.0 

6.8 

6.3 

4.1 

16.0 

.350 

9.0 

23.3 

8.6 

8.3 

6.0 

25.8 


AVERAGE OF 3 FIELDS. 15 TO 64 SPECTRA 
TAKEN AT EACH DATE 


BAND AND DATE; HAND COUNTY. S.D. 


OATS cn 


4/21 

5/10 

6/01 

6/16 

7/07 

7/27 

13.4 

9.3 

18.8 

29.1 

31.4 

11.6 

14.8 

8.1 

15.9 

29.3 

21.0 

9.0 

13.6 

7.7 

14.8 

31.8 

17.5 

10.9 

10.3 

5.6 

16.3 

31.1 

17.2 

12.2 

3.9 

4.1 

12.6 

37.0 

24.1 

5.2 

5.1 

4.0 

23.6 

40.0 

41.9 

12.0 


AVERAGES OF 2 FIELDS, 57 TO 77 SPECTRA 
TAKEN AT EACH FIELD 





BAND DATE: HAND COUNTY, S.D. 


PASTURE Cy% 


4/21 

5/10 

6/01 

6/16 

7/07 

7/27 

15.3 

15.2 

17.7 

20.8 

15.0 

13.4 

12.3 

11.7 

11.5 

15.6 

11.3 

10.1 

11.1 

n.o 

10.6 

15.1 

10.5 

9.7 

11.5 

9.2 

12.7 

21.8 

9.7 

9.3 

13.3 

10.2 

14.1 

18.7 

10.1 

10.1 

15.0 

16.1 

24.1 

26.8 

15.7 

14.9 


AVERAGES OF 4 FIELDS, 76 TO 142 SPECTRA 
TAKEN AT EACH DATE 




APPENDIX E 


INFOBMATION EXTRACTION FROM RADIOMETRIC DATA 
DISCRIMINATION THEORY 


There is a distinction between identification and dlscrlininatlon. 
DlBcrlmination means separation among species, i.e. their grouping into 
distinct classes, plus the assignment of a' numerical value to each class: 
in our case, the total ground surface occupied by that class. The func- 
tion of identification requires in addition the ’’naming** of each class, 
l.e. its assignment to predetermined categories. 

The function of discrimination requires the sufficiently accurate 
recognition of distinctions between classes: Identification requires in 

addition the Identification of specific characteristics which univocally 
tag each distinct class as pertaining to a preassigned set. 

The treatment which follows address-'s the function of discrimination 
only, specifically: 1) the identification of the drivers of accuracy} 

and 2) the statement of their interrelationships. 

Single Band - Binary Mixture Case 

Consider an area containing total pixels of two crop species: N^ 
of species A and of species B. The reflectances of species A, B are 
characterized, by Gaussian distributions with mean reflectances y., p 

A ii 

and standard deviations a., a_. 

A B 


. F,-2 


The proportion of the area containing species A is: 


o ■ N./N 
A 


That containing species fi Is; 


(1 - a) - Ng/N 


The probability density function p (x) of the distribution of the 
binary crop mixture is; 


p (x) “ <* (x) + (l-a) Pg(x) 


(E-1) 


where Is the conditional probability of x belonging to Class A and 
Pg is the conditional probability of x belonging to Class and a Is the 
a priori probability -of n belonging to Class A. 


The expected mean reflectance of the N pixels Is given by: 


X p (x) dx 


to 


X [ot (x) + (1-a) Pg (x)] dx 


ct + (1-a) pg 


(E-2) 


and the expected variance by: 


E-3 



2 

X p(x) dx “ 



7 ? [a (x) + (1-a) Pg (x)I dx - 


« OLO^ + (1-a) Og^ + a (1-a) (pg-y^)^ 


(E-3) 


For sufficiently large sasiples the sample mean M tends to a normal 
variable with mean y^ and standard deviation Using the standard 
definition of the q% probablllf/ that the true value is contained in an 
Interval 


M - k (-^) < y < M + k (-^) 


where t is a constant whose values are given below: 


CONFIDENCE LEVEL OF 


LEVEL SIGNIFICANCE 


q% 

(100-q)% 

k 

90 

10 

1.64 

95 

5 

1.96 

99 

1 

2.58 

99.9 

0.1 

3.29 

99.99 

0.01 

3.89 


For species A with y^ as the estimate of the mean value for species A, 

the true value y. lies between the limits 
A 


E-4 


“a- 


k o. k o. 


(E-4) 


/N. 


N. 


at the confidence level q. Similarly for species B: 
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Adding the two equations and weighting them by their a priori 
proportions, yields: 


aiiB + (1”“) Og - k (a + (l-o) + (1-a) Pg < 

A D 


< aii, + (1-a) P_ + k [a + (1-a) ] (E-6) 

A o 

Rewriting the equation: 


0 (1-a) o 

M - k [a 


A B 


< M + k [a 

vC *C 
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Where: 


M « a + (1-a) pg 


Introducing the definitions below In (E-7) : 


a 





a a. 


yields: 


5 [ /~a + b / 1-a ] < a 

a 


< g + _iS [ 4. jj 

av^ 


(E-8) 


Thus the absolute value of the error e of the estimate of ci 
Is bounded by: 


e < 


*¥ 


[/~oT + b /l-a ] - 


J ^ 

/ a a* + 


'‘A - “b 


A ■ °t 


a 


(E-9) 
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